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Abstract: We develop a Boltzmann- type quantum transport theory for interacting 
fermion and scalar fields including both flavour and particle-antiparticle mixing. Our 
formalism is based on the coherent quasiparticle approximation (cQPA) for the 2-point 
correlation functions, whose extended phase-space structure contains new spectral shells 
for flavour- and particle-antiparticle coherence. We derive explicit cQPA propagators and 
Feynman rules for the transport theory. In particular the nontrivial Wightman functions 
can be written as composite operators ^ AFA, which generalize the usual Kadanoff-Baym 
ansatz. Our numerical results show that particle-antiparticle coherence can strongly in- 
fluence CP- violating flavour mixing even for relatively slowly- varying backgrounds. Thus, 
unlike recently suggested, these correlations cannot be neglected when studying asymmetry 
generation due to time-varying mass transition, for example in electroweak-type baryoge- 
nesis models. Finally, we show that the cQPA coherence solutions are directly related to 
squeezed states in the more familiar operator formalism. 
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1. Introduction 



A new quantum transport formalism for coherent, interacting quantum fields based on a 
distributive expansion of the dynamical 2-point correlation functions was recently intro- 
duced in refs. [1-6]. Indeed, it was observed in [1] that in systems with certain space-time 
symmetries the information of nonlocal coherence is encoded in new spectral shell solutions 
for the dynamical 2-point functions at particular off-shell momenta. In the leading-order 
expansion in gradients this singular structure is called the coherent quasiparticle approxi- 
mation, or cQPA. It allows writing down a self-consistent network of Boltzmann-like trans- 
port equations for moments of the correlation functions, or equivalently, for the phase-space 
distribution functions. 

So far a complete formulation of the theory has been given only for non-mixing scalar 
and fermionic fields [6]. However, in nature most of the interesting applications involve mix- 
ing fields, i.e. fields whose mass-eigenstates do not coincide with interaction eigenstates. 
This certainly is the case for neutrino oscillations in the early universe [7], electroweak 
baryogenesis (EWBG) [8-15], models for spontaneous (or coherent) baryogenesis [6,16], 
and for variants of leptogenesis [17-20]. In this paper we extend our formalism to the case 
of flavour mixing. The result will be a quantum transport theory including both coher- 
ent flavour-mixing and coherent particle-antiparticle-mixing effects in temporally varying 
backgrounds. In the companion paper [21], we present a flavour-covariant formulation 
of the formalism, including also nonzero dispersive self-energy corrections giving rise to 
true quasiparticle dispersion relations. Other recent works in a similar context include 
e.g. [14,15,22,23]. 

In refs. [1,2] the fermionic transport theory was formulated using a less familiar spin- 
projected 2-component notation. Scalar fields were considered in [3] and the need for 
resummation of the collision integrals involving the coherence functions was pointed out in 
refs. [4,5]. A complete formulation of the theory in terms of the familiar 4-component Dirac 
notation and with a compact set of Feynman rules for non-mixing fields was presented in 
ref. [6]. While many aspects of the current derivation are similar to the non- mixing case, the 
added complication of the resulting equations is nontrivial. Also new physical phenomena 
emerge, such as new mixing sources in the quantum Boltzmann equations (qBE) coming 
from the gradients of the diagonalizing matrices. As in ref. [6] we will only consider the 
spatially homogeneous and isotropic problems in this paper. 

We will derive a multi-field extension of the momentum space Feynman rules presented 
in [6]. Moreover, we express these rules in a new form which is superior to the old one in 
clarity and ease of use, also in the case non-mixing fields. In particular the vertex rules 
of [6] involving complicated energy signatures are now absent. The only non-standard 
Feynman rule concerns the Wightman functions which acquire a composite form: G^^^ ~ 
AaF^i)'^Ab' Here A is the spectral function and the effective 2-point interaction vertex i^^'^ 
encodes all coherence information between the mass-shell states a and b of possibly different 
flavours and/or particle- antiparticle nature. These rules enable a systematic diagrammatic 
evaluation of the self-energies and eventually the collision terms of the Boltzmann-like 
transport equations for flavour-coherent fermionic and scalar fields, valid for arbitrary 
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types of interactions. Whenever a specific model is needed, we will use the Yukawa theory 
as an example. In the present work, the homogeneous background is always modeled by 
a time-dependent mass matrix. See ref. [21] for dispersive self-energy corrections and true 
time-dependent quasiparticle dispersion relations. 

As a numerical application we solve our cQPA induced qBE's for two mixing fermionic 
fields with a complex, time- varying mass matrix in the presence of decohering interactions. 
We study in particular the generation of the particle-antiparticle asymmetry through the 
CP-violating terms in the mass matrix. A similar mixing-scenario for scalar fields was 
recently considered in ref. [14], where the particle-antiparticle mixing solutions were ne- 
glected in the qBE's and the resummation procedure of the collision terms was not carried 
out. Surprisingly, and in contrast with the claim in ref. [14], we find that neglecting the 
particle-antiparticle correlations is not warranted despite the rapid oscillation frequency 
~ 2a; in this sector, even in the regime of relatively slowly varying mass-profile. We show 
explicitly that particle-antiparticle mixing leads to a significant contribution to the large 
momentum modes of the flavour off-diagonal correlator. These contributions persists to 
late times even though the particle-antiparticle-mixing correlations are present only for a 
short time during the transition period (when the mass matrix is varying). The integrated 
contribution from these modes is typically dominant over the asymmetry created without 
the particle-antiparticle mixing, and therefore these effects should be very important for 
quantitative calculations of EWBG-type scenarios. 

Finally, we show that the singular cQPA coherence solutions are related to squeezed 
states [24] in the operator formalism language. Indeed, by using a complex scalar field as 
an example, we first show that flavoured squeezing corresponds to flavoured Bogolyubov 
transformations (providing a generalization of the standard non-flavoured case) and then 
that there is a one-to-one correspondence between the on-shell functions flj±^ ffj± and the 
independent parameters of the flavoured Bogolyubov transformations. A similar connection 
between squeezing and coherence was recently found in ref. [22], where entropy production 
in out-of-equilibrium quantum field theory was analyzed. 

This paper is organized as follows. In section 2, we review the derivation of the 
fermionic cQPA shell structure and the equations of motion which follow from using the 
cQPA structure as an ansatz in the full dynamical Kadanoff-Baym equations. We also in- 
troduce a convenient parametrization of the flavour-coherent propagators using the familiar 
4-component Dirac notation with projection operators. In section 4, we present the resum- 
mation method giving rise to a consistent expansion of the fermionic collision integrals, 
and we give the generic forms of these integrals for spatially homogeneous and isotropic 
problems with arbitrary types of interactions. In section 5 we repeat the analysis for scalar 
fields, and in section 6 we present the generalized Feynman rules for the computation of 
the resummed self-energies for the Yukawa-theory. In section 7, we present the numerical 
example of two-fermion mixing, and in section 8 we show the connection between cQPA 
flavour mixing and squeezing in the operator formalism. Finally, section 9 contains our 
conclusions and outlook. 



-3- 



2. cQPA for fermions 

We shall first briefly review the basic theoretical setup here. We will use the Schwinger- 
Keldysh (or Closed Time Path (CTP)) approach [25] to nonequilibrium quantum field 
theory (see e.g. [26]). The main quantities of interest for us are the Wightman functions 
iS'^{u,v) = {i/j{v)ilj{u)) and iS^{u^v) = {ilj{u)'ijj{y))^ ^ which in the Wigner representation 
become: 

S{k,x) = y*dVe^^-^5(x + §,x-§), (2.1) 

where x = {u^v)/2 is the average coordinate, and k is the internal momentum conjugate to 
the relative coordinate r = u — v. These functions obey the Kadanoff-Baym equations [27] 
(see also [2,12]): 

(^ + ^ -mo - im57')5<'> - e-^0{S^}{^<'>} - e-'0{E<^>}{S^} = iCcoii , (2.2) 

where = - {S> - S<)/2 and = - (S> - S<)/2, where and denote the 
time ordered Green's function and the corresponding self-energy, and the collision term is 
given by 

Ccoii = ^e-^ ({S>}{5<} - {S<}{5>}) . (2.3) 
The (>-operator is the following generalization of the Poisson brackets: 

0{f}{9} = \ [dxf ■ dug - dkf ■ d,g] , (2.4) 
and the mass operators mo and ms are defined as: 

mo,55(fc, x) = mh,a{x)e-l^--^^Sik, x) . (2.5) 

Here rrih = {m + m))/2 and rria — {m — m))/{2i) are hermitian and antihermitian parts of 
a complex and possibly spacetime dependent N x N mass matrix m in flavour indices. The 
components of self-energy S are in general complicated functionals of the correlators 5^'^, 
which need to be computed using some truncation scheme, e.g. 2PI effective action [28] 
in loop expansion. The spectral function = ^{S^ + S^) obeys an almost identical 
equation to Eq. (2.2), but without the collision term on the RHS, and S"^'^ replaced by 
F = |(S^ + S^) in the last term on the LHS of Eq. (2.2). Thus, in the non-interacting 
limit the equations of motion for the spectral function and for the dynamical Wightman 
functions are the same. However, the spectral function is further subjected to the spectral 
sum-rule: 

J^A^ik,x)^' = l, (2.6) 

which follows for example from the equal-time anticommutation relations of the fermionic 
fields. We will see below that this condition is strong enough to completely determine the 
form of the free-theory spectral function. 

^Note that our definition of differs by sign from the more standard convention 
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2.1 Approximations 

Compared to standard kinetic approach [12,26,29] the key idea of coherent quasiparticle 
approximation (cQPA) [1-6] is to give up the assumption of nearly translation invariant 
2-point correlators. This allows new singular shell solutions, which are oscillatory in space- 
time at quantum scales ~ fc, and describe nonlocal quantum coherence between the usual 
quasiparticle excitations. In order to find cQPA we make the following approximations for 
the KB-equations (2.2): 

1) We neglect the terms oc S^. 

2) We neglect the terms oc S^. 

3) We consider only the spatially homogeneous and isotropic case. 

The first approximation is made in the standard kinetic theory as well, and it should apply 
in the limit of weak interactions when the finite (scattering) width effects can be neglected. 
The second approximation is made here only for simplicity. These terms would give rise to 
modified dispersion relations for the quasiparticles, but not change the generic structure of 
the theory. We consider the case with a nonzero dispersive self-energy in the companion 
paper [21]. The third approximation with dxS'^'^{k, x) = dxm(x) = is essential, since the 
simple spectral structures for the coherence solutions arise only in systems with particular 
space-time symmetries [Ij. With the assumptions above, the KB-equations (2.2) can be 
reduced and decomposed into the following hermitian (H) and antihermitian (AH) parts: 

(H) : 2koS< = HS< + S<H^ + i^^C^n - Cj,u)7° , (2.7) 
(AH) : idtS< = HS< - S<H^ + i7°(Ceoii + C^U)7° , (2-8) 
where we have defined a hermitian Wightman function S'^ = iS'^j^ and the operator 

H = k'd + 7^ mo + i^^^^rh^ , (2.9) 

which can be viewed as a local free-field Hamiltonian operator in the presence of time- 
varying mass. Note that only the AH-equation contains an explicit time derivative acting 
on . Thus Eq. (2.8) is mainly controlling the time evolution of the Wightman function 
5^, and it will be referred to as the '^kinetic equation^\ The H-equation (2.7) on the other 
hand, becomes algebraic in the collisionless limit and in the lowest order of time-gradients. 
The role of this ^'constraint equation!^ will be to restrict the phase space structure of . 

There are two more approximations to be made before we arrive to the cQPA spectrum 
and equations of motion: 

4) We expand the constraint equation (2.7) to the zeroth order in the scattering width 
r and in the mass gradients dtm. This gives rise to the singular cQPA phase-space 
structure. 

5) We insert the singular cQPA structure as an ansatz to the dynamical equation (2.8), 
and expand this equation to the lowest nontrivial order in F and the mass gradients.^ 

^The actual dimensionless expansion parameters are roughly dtm/uP' and V /uj. See ref. [21] for further 
discussion. 
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The set of approximations 1-5 comprises just the standard approximations in the kinetic 
theory. However, in cQPA we do not assume (approximate) translation invariance of the 
2-point functions S{k^t). We shall see in the next section how these assumptions give rise 
to the singular cQPA phase-space structure with the on-shell distribution functions fj^^^ 
and f^^^ containing SA/"^ independent real components for N x N mixing in flavour space. 
When this structure is fed into the dynamical equation (2.8) and the equation is integrated 
over fco, we find a closed set of extended quantum Boltzmann equations for the on-shell 
functions fa = f^h±' 



2.2 Phase-space shell structure 

Let us now analyze the constraint equation (2.7) in the approximation defined above, i.e. to 
the lowest order in dtm and neglecting the collision terms, which gives: 

2koS< = {H,S<}, (2.10) 

where {a, b} = ab + ba is the usual anticommutator and 

H = k'd + ^^rrih + i^^^^rria . (2.11) 

We choose to work in the mass eigenbasis, and so we first diagonalize m by a biunitary 
transformation: 

UmV^ , (2.12) 

where U and V are the unitary matrices that diagonalize the product matrices mm^ and 
m^m^ respectively, chosen such that = have real mass eigenvalues m^. With these 
definitions the direct product matrix (here (g) separates the Dirac and flavour indices) 

Y = Pl®U + Pr®V, (2.13) 

diagonalizes the Hamiltonian: 

= YHY^ = k . 5 + 7^md . (2.14) 

Multiplying Eq. (2.10) from the left by Y and from the right by Y^ we then get an identical 
equation in the mass eigenbasis: 

2A:oS'^ = {i/d, S^} , where = YS^Y"^ . (2.15) 

Next, we notice that the most general spatially homogeneous and isotropic 2-point function 
S^{k,t) can be parametrized in terms of helicity projectors: P/^ = ^(1 + /ik • 7^77^) with 
k = k/|k| and h zbl, as 

(^^ t) = lj2^h [^ho -hk-agh3 + l^Qhi - il^l^9h2] , (2.16) 
h 

where gha{k,t) are hermitian N x N matrices in flavour indices. This is a convenient 
parametrization for the problem, because the helicity operator commutes with the Hamil- 
tonian H (and with the transformation matrix y), implying that helicity is conserved in 
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a collisionless theory. Using the chiral representation of the Dirac algebra: 7^ = 1, 

7 = icr^ (g) a and 7^ = —a^ (g) 1, the hehcity block-diagonal Wightman function (2.16) 
becomes 

St = Y.3h® ^(1 + hk'S), with g< = ]^ {gho + 9h'3) , (2.17) 

h 

corresponding to the standard Bloch representation for the chiral part g^{k^t). Inserting 
this decomposition into the constraint equation (2.15) gives two (for h = ±1) homogeneous 
matrix equations for g^a ' 
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( ko /i|k| —fhij \ 

/i|k| ko —iArriij 

-fhij ko 

y iArriij ko J 



(2.18) 



where the dependence on the Bloch indices (index ordering is here defined as a,/3 = 
0,3, 1,2) is written as an explicit 4x4 matrix, and we denote fhij = (m^ + mj)/2 and 
Aruij = {rui — mj)/2. These equations may have nonzero solutions only if the determinant 
of {Bh)ij vanishes. This 4 x 4-determinant is easily evaluated for each flavour element ij, 
giving rise to N'^ independent constraints: 

det [{Bh)ij] = {kl - k^ - Mij)kl + {Aruijrhijf = , (2.19) 

where we further denote Mfj = (m? + m^)/2. The roots of Eq. (2.19) define the nontrivial 
dispersion relations: 



ko = + ^j) = =t^^J 5 or ^0 = ~ ^j) = 5 (2.20) 



where uji = (k^ + m?)^/^. These conditions give rise to a singular shell structure: g^a ^ 
6{koT^ij) or 9ha ^ ^{koTAuJij). The exact matrix form of the singular solutions can then 
be worked out by using the matrix equation (2.18). We will consider flavour diagonals and 
off-diagonals separately. 

2.2.1 Flavour diagonals 

The analysis for flavour diagonals is equivalent to the single-flavour case, considered in 
ref. [6]. Here we present it for completeness, and we also choose a different normalization 
for the mass-shell distribution functions for later convenience. For i — j the dispersion 
relation (2.20) reduces to ko = ±0;^ and a double root at ko = 0. For ko = ±0;^ Eq. (2.18) 
gives relations (we suppress the flavour indices for clarity in the intermediate steps): 

k m ' 

9h3 = -h-—gm , 9hi = -r9m , 9h2 = 0. (2.21) 
ko ko 

The remaining equation for g^o then becomes just {k^ — k^ — rrif)gho = 0, which has the 
spectral solution: 

{9ho)ii = 27rJ2±—fuh±S{koToJi) . (2.22) 
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For /co = we get the conditions 

777/ ' 

9hs = ^^^9hi , gho = 0, (2.23) 
while Qhi and gh2 are free functions at /cq = 0. This solution can then be parametrized as 
idkih = 27r^ Yl &^(^o) , {9H2h = 27r^ ^fm±S{ko) , (2.24) 

where the new subscript c refers to coherence. The reason for this particular parametriza- 
tion is that in this way the on-shell functions fi^±i}^^t) will satisfy simple oscillatory 
zeroth-order equations of motion. Using Eq. (2.16) and Eqs. (2.21-2.24) the combined 
4 X 4-spinor matrix solution is finally found to be (an identical expression holds for Sf^) 

tS<, {k,t) = 27rJ2 ±^Ph{p,^+m,) [± ^/-< 5(feo Tc.,) + (7° T ^) mjih)] , (2.25) 



where we defined fc^ = (zbcc;^,k). We have dropped the index with the understanding 
that we are always working in the mass eigenbasis when the flavour-indices i, j are explicitly 
written down. The first term in the brackets is the standard mass-shell solution with 
dispersion relation ko = icj^, while the second term represents the new ko = 0-shell 
coherence solutions. This solution corresponds to the propagator for the single flavour case 
derived in ref. [6], except here we have used a non-standard normalization for the mass- 
shell distribution functions fl^f^O^^t) for later convenience. These distribution functions 
correspond to the phase-space densities for particles and antiparticles with mass rrii and 
helicity /i, via Feynman-Stuckelberg interpretation: 

TD ■ Tn ■ 

n^^H{t) = -fa%i^,t) , n,Mt) = 1 + -7/r.< (k,t) , (2.26) 

while the new /cq = 0-shell functions /^^^^^(k, t) measure nonlocal quantum coherence be- 
tween the (flavour-diagonal) positive an negative energy modes. Finally, we observe that 
in the limit of thermal equilibrium: i^/^^^^ ^ /eq(±^z) = {e^^^' + 1)"^ and ^ 0, 
Eq. (2.25) reduces to the standard thermal correlator: 

^^zteq(^) = 27rsgn(A:o)(^ + m,) U^{k^) 6{k^ - m]) . (2.27) 
2.2.2 Flavour ofF-diagonals 

By assuming ^ rrij^ the dispersion relations (2.20) give unique solutions for equation 
(2.18), because now one can always divide by /cq ^ 0. Equation (2.21) then generalizes to: 

Qhs = -h n^i^^ ^ /^o ' 9hi = -^9ho , 9h2 = '—gho . (2.28) 

The function g^Q now gets support on all branches of the dispersion relations (2.20). It can 

be parametrized as 

{9Ho)ij =^J2m'^ + '^)&{k,^^i3) + - "^^mjiko T Aa;,,)] , (2.29) 



^We will find below that the limit of degenerated masses: rrii rrij, is well defined and smooth. 
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Figure 1: Graphical representation of the Dirac structures of different energy-components of the 
Wightman functions iSf^^ {k^t). Three-momentum and fermion number flow from left to right. 
Arrows below and above the fermion lines represent on-shell energies. Here ]Ci± = ji.^ + and 
all normalizations have been omitted. 

where the particular normalizations were chosen to facilitate the comparison with the 
diagonal limit. Inserting the results (2.28-2.29) back into Eq. (2.16) we find the full 4x4- 
spinor solution to be (again an identical expression can be written for S^j): 



h± ^^^^^ 



- + mj)f^<^5{ko T Auij)] . (2.30) 



The spinor structure of the propagator (2.30), illustrated in Fig. 1, is very suggestive: 
in addition to explicit helicity projectors it consists of projection operators -\- m 
on states with positive and negative energies for both indices i and j. It is clear that 
mixes different energy and flavour eigenstates. The on-shell distributions fl^^^ with 
i ^ j and f^Jj^^ are then naturally interpreted as describing the amount of flavour and 
particle-antiparticle coherence between the mass eigenstates with energies zbcj^ and icjj, 
respectively. In what follows, we shall use the following naming convention for the different 
contributions to the cQPA Wightman function 5^: 

"mass-shell" fl^<^, ko = ±uj, 

"particle or antiparticle flavour-coherence" fijh±^ = i^ij, ^ 7^ j 

"particle-antiparticle coherence" fiih±^ fco = 

"particle-antiparticle flavour-coherence" fij^±i = iAcj^j, i ^ j 

The complete dispersion relations corresponding to this shell structure are illustrated in 
Fig. 2 for the case of two-flavour mixing. Upper and lower blobs on the right side indi- 
cate usual flavour- mixing scenario for particles and antiparticles, including the mass-shell 
states and the particle (or antiparticle) flavour-coherence states. Middle blob corresponds 
to particle-antiparticle flavour- coherence states. In section 8 we show how the different 
coherence solutions emerge from the operator formalism. 

Finally, we note that in the limit of mass degeneracy: rui ruj^ the solution (2.30) 
actually reduces to the flavour-diagonal propagator (2.25). The limit of the propagator is 
thus correct despite the fact that the first step in the derivation is not a legitimate operation 
for ko = 0. Hence, Eq. (2.30) is the complete solution for the Wightman function Sf-{k^t) 
in the cQPA for all energies and all flavour indices. Now also the reason for choosing the 
particular normalizations is evident. Furthermore, the hermit icity of S^{k,t) implies that 
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Figure 2: Shown is the cQPA shell structure for the case of two-flavour mixing. Heavy state 
1 (blue) has mass mi = 3 (in arbitrary units) and light 2 (red) m2 = 2. Each curve is labeled 
by its energy eigenvalue and all particle-antiparticle coherence solutions are shown with dashed 
(green) lines. The number in parenthesis beside each eigenvalue gives the (real) degeneracy of the 
corresponding solution. 



the complete flavour matrices obey the hermiticity conditions 

fr±^ = fr±^ and r,t^ = r,<. (2.31) 

2.3 Spectral function and the pole propagators 

The spectral function = ^{>S^ + S'^) satisfies an identical constraint equation to 
Eq. (2.10) and thus the solution for a mass basis zj-element is identical to Eq. (2.30). 
However, the spectral function must in addition obey the sum rule (2.6), which in the mass 
basis reads: 

J ^4(fc,t)7° = <5.,. (2.32) 

This additional constraint fixes the spectral on-shell functions to fljf^ = ±^^6ij and 
fijh± ~ ^' such that the spectral function reduces to the usual vacuum expression: 

At{k,t) = ^sgn(A:o)(^ + m^)S{k^ - m])^,^ . (2.33) 

This result emphasizes the fact that quantum coherence is a dynamical phenomenon; it 
is entirely contained within dynamical correlation functions, while the static phase space 
structure of mass eigenstates defined by the spectral function, remains unchanged. This 
division is crucial also for understanding why the (flavoured) Kadanoff-Baym ansatz: 

iS<^{k,t)^2fi;:<{k,t)AUk,t), (2.34) 
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which is typicahy assumed in kinetic theory, completely fails in describing nonlocal quantum 
coherence. Indeed, the KB-ansatz explicitly forces iS^ to have the same phase-space 
structure as the spectral function, setting all coherence functions and fl^f^ to zero. 
This assumption is simply not warranted in general, because is in no way subjected to 
the spectral sum rule. However, the relation —2iA^ = + together with expressions 
(2.30) and (2.33) does give additional constraints between the distribution functions: 

-pni> \-!^^.. _ fC> _ _fC< frx qr\ 

Jijh± ~ Jijh± ' Jijh± — Jijh± ' l^Z.^Oj 

llli 

These conditions, unlike the Kadanoff-Baym ansatz, do not exclude the coherence solutions, 
but merely express Sf- in terms of Sf- . 

Let us conclude this section by noting that the retarded and advanced propagators are 
directly related to the spectral function: S'^'^{u^ v) = =f2z 9 [ ± {uq — vo)]A{u^ v)^ and so, in 
this approximation scheme, they are given by the standard vacuum expressions (e > 0): 



rrii 



However, the Feynman and anti-Feynman propagators have nonvanishing coherence parts, 
which are explicitly given by the relations 

Sjj{k^t) = Slj{k^t) — S^-{k^t) ^ 

S!j{k,t) = -Sfj{k,t)-S<^ik,t). (2.37) 

Based on Eqs. (2.33-2.37) we conclude that the on-shell functions fl^f^^ form a complete set 
of dynamical variables for the cQPA, and we will only need the kinetic equation for to 
derive a closed set of equations of motion for them. From now on we will set fl^f^^ fljh± 
whenever there is no risk of confusion. 



3. Equations of motion for fermions 

To derive the equations of motion for the on-shell functions f^j^j^^ we need to transform 
the kinetic equation (2.8) to the mass-diagonal basis. Because the mixing matrices U and 
V are in general time dependent, this procedure introduces derivative terms: 

Y{dtS<)Y^ = dtS< - , S<\ = DtS< , (3.1) 

where E'{t) = iY{t)dtY^t); we use the explicit superscript ' to indicate that is propor- 
tional to the time derivatives of mass matrix m. Using this notation we find the transformed 
kinetic equation for ij-element of the mass basis (we again drop the explicit 6/-indices): 

DtS<j = -i{H^S<^ - S<^h]) + j^Cij + 4)7° , (3.2) 

where Cij = {XCco\\X^)ij with X = j^Yj^^ and the mass eigenbasis Hamiltonian 

Hi = k-a + 7Vie-t^*^*o (3.3) 
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now involves the covariant derivative Dt in the exponent instead of the standard df. We 
now insert the spectral solution (2.30) into Eq. (3.2) and integrate over /cq, which gives 



where we denote (•••) = J ^ ( • • * )? defined (S) = S and 



-^[i^eff,<s<]^,•+7^c,,■+4)7^ 



(3.4) 



(3.5) 



Note that the fco-derivatives in Hi have now vanished as total derivatives, reducing Hi 
to the familiar diagonal form given by Eq. (2.14). The explicit form of the integrated 
Wightman function is: 

S<^=J2 PhPiil' {Pj±f^H± + PjTnjh±) ^ E i^^± + ^V^±) ' (3.6) 

h± ± 

where we have defined positive- and negative-energy projectors onto mass eigenstates: 



2 V ujj / 



(3.7) 



Using the identity Hf = uf it is easy to verify that Pi±'s are projectors obeying: 



pr± = Pi±, 



Pi+Pi- = Pi-Pi+ = and HiPi+ = ±UiPi 



(3.8) 



3.1 Flavoured quantum Boltzmann equations 

Equations (3.4-3.5) are in many ways the simplest and most compact form of the kinetic 
equations in the cQPA scheme. However, it is also useful to derive explicit equations for 
the on-shell functions. The easiest way to do this is to first solve f^f^^ from Eq. (3.6) by 
taking projections and tracing over Dirac indices: 



m 

ijh± 



nC 

J ijh:L 



iV§TV 



PjTl''Pi±PhS<j 



(3.9) 



where the normahzation constants are 

Arm ^ 

mij 



J Arc ^i^j 

and N^. = ^ 

CIJ 



(3.10) 



and we defined 



m: 



(3.11) 
(3.12) 



The equations of motion for /^^/^^ can now be obtained by taking time-derivatives of 
Eqs. (3.9), and then using the kinetic equation (3.4) for dtS^- in the trace, the projec- 
tion identities (3.8), as well as the result dtHi — 7^m^, and finally computing the traces 
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of Dirac algebra. It is convenient to define the following generalized (anti)commutator 
operations in flavour space: 

[o, z,]'" = ofs - {Ofs)^ , {o, /,r = Ofs + (0/,)t , 

[O, fs] ' = Ofs- {OUsV , {O, fsV = Ofs + (0/_,)t , (3.13) 

where O and fs are matrices in flavour indices with a (possible) dependence on a generic 
energy index 5 = ib.^ Using this notation we flnd the following coupled set of equations of 
motion for f^.^^■. 

dtfl^h± = Tt2Au,, fl-^^ - -^{^, 4-±}^ + + C;r±[/].. , (3.14) 

mij 



dtf^jH± = Ti2u,ijffjh± - -^{^' fhT^iJ + ^^h±[fh + > (3.15) 



"cij 

where the X^_^[f]-terms involving all flavour-mixing gradients are given by 



mij 

mij 



and 



([s'+,.k4-^]:,TMs'-,4V?,), 



^ 



the diagonal matrix kernels are respectively 



= '^\^ij = ^ij— and (— j^^- = d^^— , (3.18) 



and we have used a further shorthand notation: 

4\ = ^(4"± + 4%), and -4^^). (3.19) 

Finally, are the hermitian chiral components of the iEi'-matrix appearing in the com- 



mutator in Eq. (3.1): 



K^^'oiydtV^^UdtU^).., (3.20) 



I J 2 ^ 



^Note that |{0, f}^ is just the hermitian part of matrix Of and — f [O, f]^ is the anti-hermitian part 
of same matrix. The same interpretation does not apply for the [O,/]^ and {O, /}^ constructs because of 
the change of the energy index in the second term. 
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and the collision integrals are given by: 



CiAfh = iV§Iv[p,-^7°^i±A7°(Ci,- + 4)7°] . (3.21) 

Eqs. (3.14-3.15) together with the collision integrals (3.21) form a closed set of equations 
of motion for the on-shell functions /q,, once the self-energy functionals S^'^ appearing in 
the collision integrals are specified. Note that despite their somewhat involved appearance, 
they still are of a similar form as the usual kinetic Boltzmann equations, only with a larger 
number of independent variables. Also, note that most of the complications arise from the 
explicit rotation to the mass eigenbasis; the starting equation (3.4) is remarkable simple in 
the matrix form. More explicit expressions for the collision integrals will be derived in the 
section 4.2 below (see equations (4.17-4.18) and (4.21-4.22)). 

3.2 Single-flavour and flavour-diagonal limits 

It is instructive to consider the single-flavour and diagonal limits of the full evolution 
equations (3.14-3.15). It is easy to see that if we drop the mixing gradients and flavour 
mixing by the collision term, diagonal equations with i — j reduce to the single-flavour 
equations presented in ref. [6] with a positive real mass m = ttir. However, for a more 
complete comparison we have to account for the complex mass m = |m|e^^ in the single-field 
case. In the current approach the complex phase has to be removed by a chiral rotation 
Ua = exp(— i^7^^). After applying this rotation to both constraint and kinetic equations 
we get the following rotated single-flavour quantum Boltzmann equations: 

dtfmh± = ±l{Kh+fch+ + Kh-fch-) +Cmh±[f], (3.22) 

dtfch± = T«(2a;k + vi,he')fch± - ^KhT(fmh+ - fmh-) + VkhCch±[f] , (3.23) 

where = |k|/a;k and all dependence on the mass gradients in mixing terms is encoded 
in 

$U±^^k/^^T^^^^ (3.24) 
The new coherence functions fch± ^re related to the ones used in ref. [6] as follows: 

fch± = Vkh cos 9^fch± T sin 9^fch± + i {v^h cos 9^fch± ± sin 9^fch±) , (3.25) 

Note that equations (3.23) are even somewhat simpler than those introduced in ref. [6]. 
Beyond simplicity, they are also convenient because fch±^ turn out to be the functions 
with the (local) canonical normalization described in [6]. 

Let us now move on to the diagonal limit of the mixing case. The above example 
underlines the significance of the diagonal chiral rotations, which are relevant also in the 
mixing case. Indeed, it is important to note that the rotation matrices U and V are 
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elements of the full U{N) group. We can write the diagonalizing matrix y as a product of 
chiral C/(l)-rotations and the elements of the SU{N) groups: 



Y = Pl®U + Pr®V = Pl® U^UvUa + Pr® U'^UyUl , 



tNi 



(3.26) 



where Uy — exp(i^0y) and Ua — exp(i^^^) are the vector and axial vector U{1) trans- 
formations, and the SU{N)l^ji transformation matrices can be parametrized as U^j^ = 
exp(iaj^r^). In this parametrization the mixing-gradient matrices (3.20) become: 



+ 



1 - 

2 . 



7 



(3.27) 



The pure C/y-part of is proportional to identity operator and therefore it does not 
contribute to kinetic equations. This is as expected, since Uy is actually a pure gauge 
transformation. In single-flavour limit U^^^ = 1 and only axial vector transformation Ua 
contributes to mixing gradients in this limit, as shown above. 

It is interesting to observe that an effective (varying) complex phase is generated in the 
non-hermitian flavour-mixing scenario by the flavour-diagonal mixing gradients. Indeed, 
one can show after some algebra that the diagonal mixing-gradient terms (3.16-3.17) reduce 
to 



i 2SL J— ^ {fiih-\- ~ fiih-) + • • • 5 



2uJi 



f- 



lh± 



^ ^ Ikf 2a^ ('^^'^+ + fiih-) + • • 



(3.28) 
(3.29) 



respectively, where we have written down only the flavour- diagonal self-coupling terms. 
By identifying: — 2S^^ = 0^ + (q/^ — a'^)r!^- ^ 9\ and accounting for the normalization 
differences: v\^hf^- = fc and ±^f-^_^ = fm±, we see that the terms in Eqs. (3.28 - 3.29) 
reduce to the nondiagonal (in energy) mixing terms appearing in the single flavour equations 
(3.22 - 3.23). Since scalar-like mixing gradients do not contribute to flavour diagonals, 
this means that neglecting all off-diagonal flavour-mixing terms in Eqs. (3.14 - 3.15) reduces 
the system to just N copies of single-flavour equations, as expected. Note that taking the 
diagonal flavour-mixing gradients in consideration is necessary to get independent mass- 
function phase changes for all diagonal components. 



4. Resummed fermion collision term 



The basic quantity appearing in the equations of motion (3.14-3.15) is the zeroth moment- 
integral of the mass-basis collision term (Cd) = (XCcoU-^^), or more precisely, its projections 
defined in Eq. (3.21). In this section we present a general method of expanding this quantity 
to the first (leading) order in the interaction width and the gradient of the mass function, 
denoted by = 0{T^dtm). We follow closely the analysis presented in ref. [6] for the 
single-flavour case. 
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As was observed in ref. [6], the rapidly oscillating coherence functions fa give rise to 
a leading order term in each order of the naive (0-)gi"adient expansion of the collision 
integral: 

(Cd) = / ^X(t)ie-0({5]>(^,i)}{5<(fc,i)}-{S<(A.,i)}{5>(fc,t)})Xt(i). (4.1) 

If the coherence contributions are present only in the external propagators S^'^ , the oscil- 
latory gradients can be directly resummed in the Wigner representation [6]. In the general 
case the 0-expansion becomes very complicated however, and it is much more convenient 
to rewrite the collision term in the two-time representation: 

(Cd) = ^J dwoX{t){j:>{t,wo,k)S<{wo,t,k)-J:<{t,wo,k)S>{wo,t^^^^ (4.2) 

where the propagators are defined as 

= e-^^^^-^-<)S{K (4.3) 

J 27r 2 

and similarly for self-energies. We emphasize that we employ the two-time representation 
only temporarily to perform the resummation of the oscillatory gradients. At the end 
we will recover a Boltzmannian-type collision term involving momentum loop-integrals 
and the distribution functions /^'^(k,/;) in the Wigner representation. Next, we want 
to express the flavour-basis propagators S^'^ and the self-energies S^'-^ in Eq. (4.2) in 
terms of the mass-basis propagators, in order to use the spectral solutions of section 2.2. 
According to the definition in section 2.2, the propagators transform locally in time in the 
(fc, t)-representation, and we now define a similar transformation law for the self energies: 

S^iht) = Y{t)S{k,t)X^t), 

Sd(A;,t) = X(t)S(fc,f)yt(t). (4.4) 

It then follows that in two-time representation the transformation laws are still local in the 
average time-coordinate, and not following different end-point times wq and Wq separately: 

Sa{wo,w'„k) = r(^^^)5(«;o,u'^,k)xt(^^^) , (4.5) 
with a similar expression for T,^'^ (wo^WQ^'k). Then the collision term (4.2) becomes 
(Cd) = \l dwoX{t)X^C-^){E>{t,wo,k)S<{wo,t,k) 

-E<(i, wo, k)S>iwo, t, k))X(^±^)xt(i) 

= ^J dwo{^^it,wo,'k)S<iwo,t,k)-J:<it,wo,'k)S>iwo,t,k))+0\ (4.6) 
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where = 0{r'^ .Tdtm.rdtX).^ In the last row we have used the zeroth order term of 
the mixing matrix X expanded in mixing gradients: 

xC-^) = X{t) + dtX{t)C^) + Oid-'.X) . (4.7) 

We can now see that neglecting corrections in Eq. (4.6) is a good approximation if the 
gradients associated with the varying background, ^tlnm ^ dtX ^ and energy gradi- 
ents InS ^ uj~^ parametrically obey tcj > 1. This is consistent with our approximation 
that all the elements of the original flavour-basis mass matrix are slowly varying in time. 

4.1 Expansion of the propagators S^'^ 

Next, we want to write the collision term (Cd) in terms of matrix functions 5^^'^^'^(k, t) 
defined in Eq. (3.6), locally in the external time- variable t. By writing the Wightman 
function (2.30) as (we suppress the flavour indices ij for the remaining of this section for 
convenience) 

S< {k,t) = 2TTj2 [<5± < (k, mko T<^)+ Sp (k, t)S{ko T Aa;)] 

± 

= ^(5r(^,t)+^S:<(^,i)), (4.8) 

± 

where and 5^5 can be read off from Eqs. (2.30) or (3.6), and inserting it to the 

definition (4.3) we get the two-time representation of the cQPA- Wightman functions as 
(t=(no + ^o)/2): 

'^(^o,^o,k) = [e^^^(^o-^o)5^<'>(k,t) + e^^^^(^o-^o)5f '>(k,t)] . (4.9) 

We obviously need to expand the functions S^'^{wo, t, k) appearing explicitly in Eq. (4.6). 
However, the perturbative expressions for the self-energy functions S^'^(t, wq^ k) in general 
involve further integrations over internal vertices wq and Wq, and so we need to expand 
correlators S'd('^0 5 '^O' ■'^) with arbitrary time-coordinates wq and Wq with respect to the 
external time t — (uq + vo)/2 to zeroth order. This can be done by a simple Taylor 
expansion: 

5<,> ,)"a.»s->(*. <) . (4.10) 

The time-derivatives d^S^(k^ t) can be found from the equations of motion (3.14-3.15). To 
zeroth order these reduce to dtfj^j^ = =Fi2Aa;/^ and = =Fi2a)/^_j_, indicating that the 

on-shell functions /j^ and are oscillating in quantum scales with frequencies 2Aa; and 

^The second order correction to the equation (4.6) from the exphcit mixing gradients can be ex- 
pressed using the notation of the section 4.1 as 5{C) = | J^-^ [7°^^° , (^A;o^^(^o, k, ^)<S^(k, ^) — 
dko^ff^{ko,'k,t)S^{'k,t))]\kQ=±ujk + ••• 5 where other second order terms arise from the equations of mo- 
tion of the correlators S^'^ or exphcitly expanding in gradients of dtm,dtX. 
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20;, respectively. This holds true in the distribution sense also for the functions S'^^'^(k, t) 
and S^'^Ck^t) and we get: 

Qngm<,> ^ (^i2Aa;)^^^<'> + , 9,^5^'^ = (Ti2cD)^^f + . (4.11) 

Using these equations recursively in the Taylor expansions, we find that distributions S'^^'^ 
and S^^'^ depending on the average time (wq + Wq)/2 reduce to quantities evaluated at 
time t, multiplied by the phase factors: 

n=0 ^' 

Inserting these results to Eq. (4.3) one eventually finds: 

_^ ^■=fiAuj{wo-w[^)^iO{wo+w'Q-2t)^^,> ^-^^^^^^ _^ q1 ^ (4.13) 

As a trivial consistency check we see that this expansion reduces to Eq. (4.9) as W() uq 
and Wq ^ vq. Physically the propagator (4.13) is an approximation which takes into 
account the rapid temporal variations due to oscillations of the coherence functions, but 
neglects the corrections due to temporal variations in the background fields as well as the 
corrections from the full memory effects of the dynamical interactions. When applied to 
the case Wq = t, Eq. (4.13) implies a particularly simple phase dependence: 

S^'^{wo,t,k) = ^e^^^^(^o-^)5^'>(k,t) + , (4.14) 

± 

where we denote 5^'^ = S^^'^ + S^'^ . Applying equation (4.14) to the propagators in 
the collision term (4.6), we find 

(^^^•) = (^eV(^^^'k,t)5<^(k,t) - S<ff^,,(±c.,,k,t)5>^(k,t)) +0^ (4.15) 

k± 

where the effective self-energy function is defined by: 

Kiijik,t) ^ I d^oe^'°(*-"°)sS(i,^o,k) . (4.16) 

This is formally the most important result of this section. Note that the effective self-energy 
functions S^^^(iba;j, k, t) are not in general equivalent to the projections S^'^(d=a;j, k, t) 
of the original Wigner representation self-energy. However, both reduce to the same expres- 
sion m the limit dt^<^>(k,t) = 0\ i.e. in the case that the oscillatory coherence solutions 
do not contribute to the self-energy. Perturbation theory rules for evaluating ^^^^ij{k^t) 
will be derived in section 6. 
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It is interesting to see that in the end, after the resummation has been carried out, 
Eq. (4.15) can be written in an integral form formally equivalent to the naive lowest order 
gradient expansion of the equation (4.1): 

idj) ^Ij^ {^es,ikik, t)S^s,in,kjik, t) - S< t)5> t)) . (4.17) 

The crucial difference is the replacement of the original self-energy functions and the cQPA 
Wightman propagators by the effective self-energy functions and the effective in- 

propagators: 



c<,> 



(4.18) 



We conclude this section by finding the most general explicit expressions for the collision 
integrals Eq. (3.21), consistent with the symmetries (spatial homogeneity and isotropy) of 
the problem. 

4.2 Explicit collision integrals 

By using covariant formulation for the Dirac structures the collision integrals (3.21) can 
be rewritten as 



mij 



CIJ 



(4.19) 



Similarly to Eq. (2.16) we can write the helicity projections of the most general spatially 
homogeneous and isotropic effective self-energy in the Dirac notation as 

P,iS</ (fc, t) = Pn [7° A<'> + k • 7 B<'> + C<'> + ihj' D<'>] . (4.20) 

Now, using Eq. (4.15) for {dj) and the projector properties of the propagator matrices 

qm,c< 



^m^<,> ^j^^ self-energies (4.20) it is relatively straightforward to show that the collision 



integrals of the Boltzmann equations (3.14-3.15) take the following form: 

C±[/].. ± {Bf:,v^ ± {C>, ^ /<,-}- 

{i [A>,v^ ± i [B>,f<^]- + {D>, ^ f<^}^) 



mij 



|k|mf 



(4.21) 
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+ ( - 4<-}?,- ± 41?];. + i - 4%-]:, 
- ( [Bi , ^ 4%1 - [ci , 4V] ± ^ > 41?] ;) 



+ 



rriij 
\^\fhij 



(4.22) 



Keeping in mind that here i^^^^ij^ = depends on the energy index ib, we 

have enlarged the notation (3.13) to accommodate this: 

[5^55 fs\ = ^sfs ~ i^sfs)^ 5 {^S1 fs}^ = ^sfs + {^sfs)^ 5 

[S,, /,] ^ = S,/, - (S_,/_,)t , {S„ /J^ = S,/, + (S_,/_,)t , (4.23) 

and the generahzed anticommutators of /^:s with the self-energies S without d= indices 
are reduced as follows (with obvious generalization to commutators): 

{E,4<±'r = ^{s±,4Tr ±' li^^Jn^r, 

{s> 4±^'F = f^^r ±' ^{St, • (4.24) 

For example the first term of Eq. (4.21) is explicitly written as follows: 

k 

+ A^hi^cv,)r,<,^ + n^l^AlJ^i^cvk)) , (4.25) 
and respectively the first term of Eq. (4.22) is 

k 

+ A>./.(^^^)41^T + &All,i±oj,)) . (4.26) 

Note that the flavour index of the matrix multiplication dictates the energy dispersion 
ko = ±ujk of the effective self-energy, and observe in particular the delicate inversion of 
signs in the dispersion for the second /'^-terms, in accordance with Eqs. (4.15) and (2.31). 

Equations (4.21-4.22) present the most general form of the fermionic collision integrals 
in the cQPA approximation, consistent with the spatial homogeneity and isotropy. Fur- 
ther reduction is possible only after an explicit expression for the self-energy functions is 
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specified.^ A perturbative expansion for S^^'s, which automaticahy accounts for the re- 
summation of the oscihations of the coherence propagators, could easily be written down in 
the two-time representation using the zeroth order expanded propagators Eq. (4.13). How- 
ever, in section 6 we will show how this expansion can be rewritten in terms of effective 
momentum space Feynman rules, which make the computation of straightforward 
up to a desired order in perturbation expansion including an arbitrary number of coherent 
propagators, without leaving the Wigner representation. Before this, we shall extend our 
analysis to the case of scalar fields, however. 

5. Scalar fields 

The cQPA approximation can be employed also for other types of quantum fields than 
fermions, although the analysis is qualitatively different since the spinor structure and 
corresponding projection operators are not present. In this section we formulate the cQPA 
quantum transport equations for flavour-mixing complex scalar fields following the same 
approximation scheme we used for fermions, described in section 2.1. As a result, we will 
find the same singular phase space structure with the shells: ko = ±a; = ±{LJi ^ uj) /2 and 
ko = zbAcj = :iz{(jji — CL;j)/2, as for fermions.^ However, because of the structural difference 
between scalar and fermionic KB-equations, the derivation of the transport equations is 
now conceptually different to the method for fermions; analogously to our earlier works 
for single flavour scalar field [3,6], we need to introduce multiple fco-moinent integrals to 
obtain closure. 

As with fermions, we start by writing down the Kadanoff-Baym equations for the 
Wightman functions iA^{u,v) = {(/)\v)(/){u)) and iA^{u^ v) = (0(i/)0t(^)) in the Wi gner 
representation (see e.g. ref. [12]): 



where = - (A> + A<)/2 and H^ = H^ - (n> + n<)/2, while A^ and H^ denote 
the time-ordered Green's function and the corresponding self-energy. The form of the 
collision term C^^jj can be obtained from Eq. (2.3) by replacing S ^ A and S ^ H. The 
mass m? as well as iA^'^(/c,x) and iH^''^(/c,x) are hermitian N x N matrices in flavour 
indices. Proceeding through the steps 1-3 in the approximation scheme described in section 
2.1, we drop the terms proportional to A^ and H^ in KB-equations (5.1) and (by spatial 
homogeneity) set dxiA^'^{k^x) = dxm?{x) = 0. After these approximations we break the 

®Let us again emphasize that if there is no coherence contributions inside the self energies, e.g. in the 
case of interaction with thermal bath, then to leading order i;^^(A;,t) just reduce to standard Wigner 
representation self-energies E^'^(/c,t), which can be computed using the standard techniques of CTP- 
formalism or real-time formalism of thermal field theory. This type of simple example is considered in 
section 7, where we apply our formalism to CP- violating fiavour mixing in the presence of collisions with 
thermal background. 

^This complete shell structure for mixing scalar fields was also observed in ref. [14]. However, in this 
work the rapidly oscillating ko = zbAcj coherence-solutions were subsequently neglected in the derivation 
of the transport equations. 




e-0{nn{A<'>}-e-0{n<'>}{An = C^,, , (5.1) 
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KB-equation into hermitian (H) and antihermitian (AH) parts, finding eventually: 

(fc^ - \dt - cos (i^r^fj \ {m^ o} + sin (i^r^fj \ K, o] ) iA<'> = -Ca , (5.2) 
kodt + cos (^arafj \ [m\ o] + sin \ {m\ o}) iA<'> = Ch , (5.3) 

where [X, o\Y = [X, Y], Ch = (C^n + Cfii)/2 and = (C^^^u - cSi)/(2i), and the super- 
scripts on derivatives denote explicitly that t-derivatives always act on the mass function 
m^(t), while /cQ-derivatives act on the correlators zA^'^(/c,t). 

5.1 Phase space structure 

We now proceed to steps 4-5 in the cQPA approximation scheme, and analyze the KB- 
equations (5.2-5.3) in the zeroth order in F and dtm to find out the singular phase-space 
structure. Unlike in the case with fermions, for scalars both H- and AH-equations contain 
explicit time-derivatives even in zeroth order in mass gradients, and we need to use both 
equations to derive the algebraic constraint: first the collisionless AH-equation (5.3) is 
solved to the lowest order to give dtiA^^^ = —^[mP^iA^^^]^ and then this relation is 
used recursively to compute the second derivative term dfiA^'^ in Eq. (5.2), again to 
the lowest order in mass gradients. This procedure leads to the zeroth-order algebraic 
constraint equation 

(^o-k'-^{m2,o} + ^ [m2,[m2,o]])iA<'>(A;,i) = 0. (5.4) 

Next, we go to the mass basis by diagonalizing the hermitian mass matrix: m? = 
Um?U^ , where U is a, unitary mixing matrix. The correlators transform similarly: A -> 
Ad = C/AC/t. In the mass basis Eq. (5.4) becomes 

(^{kl - - Mf^)kl + (Am,,m,,)2)iA<.'>(A;,t) =0. (5.5) 

where Mfj = (m? -|- m^)/2. We see that the multiplying factor in parenthesis is equal 
to the fermionic constraint determinant in Eq. (2.19), and thus the nonzero solutions of 
A^'^ = (A^'^)^j must be proportional to either 6{kQ ^ Cuij) or 6{ko =F Acuij), as before. 
There are no further algebraic constraints and we can parametrize the cQPA-propagators 
A<'> as 

iA<'>ik, t) = ^J2i^ fS'^^ T Oij) + f^<^>6 (ko T Acoij) ) . (5.6) 

To zeroth order the spectral function = |(A^ — A^) obeys identical KB-equations 
to Eqs. (5.2-5.3), and consequently the solution is of the same form as Eq. (5.6). In addition 
however, the spectral function must obey the sum rule 



/ 



— (^o + k)At(fc,i) = ^^i, (5.7) 
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following from the equal-time commutation relations of the scalar fields 0^. It turns out 
that this relation completely fixes the values of the corresponding on-shell functions (see 
ref. [3]) to flj^ = ^6ij and f[j^ = 0, reducing the spectral function to its standard form 



Atj{k,t) = 7rsgn{ko)6{k^ - m^)Sij 



(5.8) 



Furthermore, using this result with the defining relation 2iA^ = — A^, one can easily 
show that the dynamic functions /^(k, t) and /^(k, t) are related: 



(5.9) 



That is, only half of the on-shell functions appearing in A^'^ are free variables. In what 
follows, we derive equations of motion for the on-shell functions //^jf^ = = fa- 

5.2 Equations of motion 

We derive the transport equations for the on-shell functions fa by inserting the singular 
cQPA-propagator (5.6) as an ansatz into the full KB-equations (5.2-5.3). The unknown 
on-shell functions can be extracted from /\f-{k^t) by integration with different weight 
functions. To be specific, we use moment functions: 



dko 
2tt 



(5.10) 



To get the closure we need at least four different moments, or alternatively, we could use 
three moments and one time-derivative of a moment since the constraint equation (5.2) is a 
second order differential equation w.r.t. t. It appears that the most natural choice is to use 
the three lowest moments po,i,2 and the derivative dtpo. In order to relate these moments, 
and in particular the derivative dtpo, to the /-functions of Eq. (5.6), one needs to solve 
the derivatives dtfa in terms of fa to zeroth order. The required zeroth order equations 
of motion can be obtained by solving Eq. (5.3) for any four moments, for example pi,2,3,45 
using the singular cQPA-propagator (5.6). In this way one finds the equations 



(5.11) 



which are identical to fermionic equations of motion (3.14-3.15) to zeroth order. By com- 
puting the moments (5.10) of the cQPA-propagator (5.6), and using Eqs. (5.11) for the 
derivatives dtfa appearing in dtPo, we then obtain a set of invertible relations between fa 
and the moment functions: 



/ PO \ 

dtPo 
Pi 
\ J 



2ujiUJj 
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-2iAa; 
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Alo -Aoj 



-£1)2 Aa;2 Atj2 J 



( f+\ 

frn 

n 

\f-J 



(5.12) 



where an element wise rule of matrix multiplication is applied to the flavour indices which 
are suppressed to generic subscripts for convenience. 
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The evolution equations for the moments po,i,2 and dtpo are found by taking the zeroth 
moment of Eq. (5.2) and the zeroth and the first moments of Eq. (5.3). Transforming these 
equations to the mass-diagonal basis amounts to replacement: 

dtPn ^ DtPn = dtPn " ^[6' , p„] , 9' = iUOtU^ , (5.13) 

analogous to fermions (see Eqs. (3.1-3.2)). The resulting equations of motions for the 
moments are given by 

D^po - 4p2 + 4(0)^ + Auj^)po = 4(CAd) , 
Dtpi + 2iAujOpo = {Cnd) , 
DtP2 + 2iAujOpi - {OO' + AujAuj')po = {k^Cud) , (5.14) 

where C(//^)(i = UCh^aU^ i and we have suppressed the flavour indices, for which the 
multiplication rule is elementwise as in Eq. (5.12). 

Unlike in the case of fermions, we do not attempt to rewrite the evolution equations in 
terms of /-functions, although it could be done e.g. by differentiating the inverted equations 
(5.12) and using the evolution equation (5.14). This change of variables carries a subtle 
issue related to the gradient expansion however, as the moment equations are exact in the 
sense of gradient expansion but the inversion relations are only valid to zeroth order. To 
avoid this unnecessary loss of accuracy, and also to keep the evolution equations simpler, 
we prefer to write the qBE's in terms of the moments rather than the /-functions in the 
scalar case. In practice, this means that the /-functions are only used in the collision 
integrals, and only in this context their relation to moments needs to be defined through 
the inverse of equations (5.12). 

The subtlety with the inversion is related to an inherent "ambiguity" in choosing the 
particular set of moments (or choice of moment weight functions) to define the integrated 
equations of motion. Indeed, transport equations based on higher moments {e.g. for pi,2,3,4) 
would be equivalent to Eqs. (5.14) only up to zeroth order in gradients and scattering 
width, because the singular correlator (5.6) and consequently the relations (5.12) between 
/-functions and the moments are valid only to zeroth order. Our choice of moments is the 
most natural one in the following sense however: using the flat weight (lowest moment) 
to integrate the KB-equations (5.2-5.3) corresponds to encoding no exterior information 
on the dynamical equations (see the discussion in section 5 in ref. [1]). This creates the 
two first equations of Eqs. (5.14). However, the fiat weight integration generates an extra 
moment p2 into the integrated H-equation, and the only way to reach closure without 
adding new moments is to take the first moment of AH-equation. This we have done to 
obtain the last equation in Eqs. (5.14).^ 

The moment choice problem is also related to which order in gradients one is work- 
ing. For example, going to higher order in gradients would involve more independent mo- 
ments, but a consistent treatment would then require using more complicated phase-space 

^For fermions this issue did not arise because there the flow term of the (/co-integrated) dynamical KB- 
equation does not couple zeroth moment of the correlator {S^{k,t)) to higher moments. For scalars the 
issue is more delicate because we are forced to use a mixture of different moment functions. 
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structure for the correlator A^(A:,t), which can be obtained as a generahzed distributive 
expansion in gradients. We explore this issue to some detail in Appendix A. 

5.3 Resummed scalar collision term 

We need to express the collision integrals appearing in Eq. (5.14) in terms of the distribu- 
tion functions and (and eventually in terms of the moments using the connection 
Eqs. (5.12)). Likewise with fermions, the basic quantity we encounter is the a-moment 
integral of the collision term in the mass basis: 

{Ca) = I ^feo"^(*)^e-^({n>(fe,i)}{A<(fe,i)}-{n<(fc,t)}{A>(fc,t)})t/t(t), (5.15) 

where we now need a = 0, 1. Again, to resum the oscillatory gradients of the distribution 
functions in the <0-expansion, we write the collision term in the two-time representation: 



{Ca) = ^ J dwo{idro) 



r-o=0 



(5.16) 



where we have dropped the time-gradients of the mixing matrix C/(t), as before. For a = 
this immediately reduces to a formula analogous to Eq. (4.6) for fermions. For a = 1 the 
extra S^-q -derivative gives rise to an additional complication, but we still continue to search 
for a consistent expansion around the external time t as before. We begin by writing the 
spectral propagator (5.6) in two-time representation as 



2uJiUJj 



where we Taylor-expanded fai}^^ wq^wq ^ around the external time variable t, and used 
the recursive zeroth order equations of motion: d^f^ = {^i2A(jj)^f^ + and 5^/^ = 
(=Fi2cj)^/^ + 0\ exactly as in the fermionic case. Now, using the expanded propagator 
(5.17) in the collision integrals (5.16) we get: 



{Ca)ij — ^ 



2 ^ 2ijJ].ijJj 



(5.18) 



±{±0,j + -a,)"in>ff .,(±a;,,k,t)/-<(k, 

+ ( ± AuJkj + ^a,)^n>,^,,(±a;,, k, t) r,<^{K t)^ 
where the time-derivative operator acts only on the effective self-energies defined as 

Kiijik,t) = J dti;oe^'=°(*-"'°)n<.'>(i,^«o,k) . (5.19) 
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Note again the simphcity of these results; to get the colhsion integral for an arbitrary 
moment equation, we only need to evaluate a single generic function n^'^^(zbct;fc, k, t), 
which is of the same form as the fermionic effective self-energy function (4.16). 

The collision integrals appearing in the equations of motion (5.14) are just the her- 
mitian and the antihermitian parts of Eqs. (5.18): {Cud) = + (C'o))' {^Ad) = 

^((Co) - {Cl)) and {k^Cud) = \{{Ci) + {C\)). Using the conventions (4.23) for the 
generalized (anti)commutators in flavour space, we finally obtain the following explicit 
expressions: 



{CAd)ij = 



± 



T ^ 



OJiUJj 



{in 



f<+i 

effdi'/di J 



<+\m 



(5.20) 



<+\m 



UJjUJ. 



(5.21) 



and 



" - ({in>ri,a;/<+}^ + u;,a;,.{m> 1 /<-}i^ ± M(m> i), /<+}^) 



+ 



+ 



(5.22) 



where the diagonal matrix kernels are defined in analogy to Eq. (3.18): 



= bi^LOi and (-)ij 
and we have again used shorthand notations: 



^12 



LUi 



(5.23) 



/<+^i(/r±/^<), /r^^(/rT/^<) and ng.^^n<^>(±a;,-). (5.24) 

As emphasized before, in order to perform practical calculations using Eqs. (5.14) with the 
moment integrals as dynamical variables, one has to use the inverse relations of Eqs. (5.12) 
to express the collision integrals (5.20-5.22) in terms of the moment functions po,i,2 and 
dtPQ. In practice, this inversion can be left to numerical routines, since the expressions of 
the collision terms take the simplest form in terms of /-functions. 

6. Momentum space Feynman rules 

In this section we derive generalized Feynman rules for computing the effective self-energy 
functions iHeff and illeff through perturbative techniques, including the coherence effects. 
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Standard methods, such as the 2PI formahsm [28], exist for diagrammatic expansion of 
the two-time self-energies T>^^{u^v) and Il^^{u^v) appearing in equations (4.16) and (5.19), 
and our task is to reduce the computation of the diagrams into a set of momentum space 
Feynman rules. We derive these rules by replacing all propagators in a given diagram by 
our resummed propagators (4.13) and (5.17), transforming these propagators to momentum 
space and performing all time integrations related to internal vertices. To be specific we 
shall consider a generic Yukawa type interaction Lagrangian between the scalar and fermion 
fields: 

Ant = -yij (f^i i^i + h.c. . (6.1) 

However, from our final results it is obvious how these Feynman rules generalize to arbitrary 
interactions. 

The main part of the derivation consists of completing the time-integrals at vertices, 
accounting for the nontrivial phase factors in propagators. To begin with, it is convenient 
to rewrite the propagators (4.13) and (5.17) in a generic 4-dimensional representation: 

A<'>{wo,w'o,P) = I ^e-^"'(-°-'='-o)+ipo(i-c.O*A<>„,^(f,,t) . (6.2) 

where S^^'^^-j{k,t) are the (effective) in-propagators, given by Eq. (4.18), with a similar 
definition for the scalar propagators ^^^n /cK^' ^)'^ phase factors are given by Cij = 
+ujj/uji for G^^^-j^ and Cij = —Uj/ui for G^^'-"^^^-, where we use generic notation G = 
{S', A} for fermion and scalar propagators. Because the phase factors are different for m- 
and c-parts of the full G^'^^^^-, Eq. (6.2) is understood to hold for the corresponding parts 
separately. For the flavour- diagonal pole-propagators we use the standard expressions, 
Eq. (2.36) in two-time representation and similarly for scalars, which are directly of the 
form (6.2) with "normal" phase factors ca — +1. We note that the normal phase factor also 
applies to the flavour-diagonal mass-shell propagators G^^'^ . For all other (component) 
propagators, which encode the flavour or particle-antiparticle coherence, the overall phase 
factors are "abnormal", except in the particular case of Wq = where the nontrivial Cij- 
terms cancel. 

To proceed, we need to show that the total t-phase arising from exp(iA:o(l— Qj)t)-factors 
in the decompositions (6.2) vanishes for an arbitrary self-energy I]eff(fc,t) or nefr(/c,t) di- 
agram. This can be shown by a direct computation using Eqs. (6.2) and following the 
identical steps as in [6] for single-flavour case. Alternatively, by the flavour-covariant for- 
mulation of [21] the cancellation of the total t-phase follows trivially by the shift of vertex 
time-integration variables. After the cancellation of the t-phase, the derivation of the 
Feynman rules is straightforward. Based on Eq. (6.2) one would obtain rules with (ef- 
fective) propagators iS^ff[^ij{k^t) and ^A^'^^ ^^(^, t) for fermion and scalar lines, and the 
vertex rule with modified energy-conservation delta-functions involving Qj-factors. These 
rules would provide a generalization of the Feynman rules derived in ref. [6] to the case of 

^In this section we use k, I instead of i,j to denote scalar flavour-indices. 
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flavour- mixing fields. However, these Feynman rules are somewhat complicated to use in 
practice, because of the nontrivial Qj-factors in the vertex rule with a delicate dependence 
on the flavour and particle-antiparticle coherence components and the direction signatures 
associated with the propagators. We shall therefore follow a different approach here. 

6.1 Extended KadanofF-Baym ansatz 

An equivalent, but more intuitive and transparent set of Feynman rules can be obtained 
by rewriting the two-time Wightman functions (4.13) and (5.17) in the doubly Fourier 
transformed form (this step is not necessary for the pole propagators): 

A<'>(^.o,^^,p) = / ^|^e-(---Pi^i)A<.> (6.3) 

where we have dropped the t-phase factor which cancels at the end, as discussed above. 
The idea here is that after performing time integrations in internal vertices wq and Wq we 
are left with normal energy conservation delta functions in momentum space at the price 
of having different 4-momenta k and at the "in" and "out" ends of each Wightman 
function. This structure is nicely reflected in the explicit forms of the doubly transformed 
propagators: 

iS<4>.{k,k';t) = Af,{k)F<l>ik,k'-t)At{k'), 

iA<il,{p,p'; t) = AUp) F<lf{p,p'; t) At{p') , (6.4) 
where the spectral functions Af^{k) and Af.j^{p) are defined in Eqs. (2.33) and (5.8), and 
F<l>(k, k'; t) ^Ai27rmk - k') ^ P,(k)4 (47^;^ + , 

h,± 

KkfiP^P'-^ i) ^<2nmp - p') ±2Co,iel (<47<'> T , (6.5) 



where e.g. 9^ = 9{±ko) and fl^^J^ = fljhJ^ 0^^^^ "^^^ ^^^zj"^™^^^^^^ understood 
as effective 2-point vertex functions, which encode all quantum statistical information of the 
Wightman functions, including flavour and particle-antiparticle coherence, as well as the 
helicity structure for fermions. Each Wightman function can then be viewed as a composite 
operator consisting of two standard on-shell spectral functions coupled by effective 2-point 
interactions; the amount of coherence in the system is directly manifested as the strength 
of these couplings. 

Now, using the following properties of the spectral functions: 



/ 



W [±^44'5i.(2vr)^5^(k-k')]4(A;') = lAt,ik)ei 



(27r)4- 

Sp' 



|4- 



Alip) ± 2u,9ielSki{27rf5Hp - p') A^ip') = Ul{p)9i , (6.6) 



2 
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we can connect our fully diagonal (both in flavour and in energy) Wightman functions to 
their usual kinetic theory counterparts: 

^ ^ ^ h,± ^ 

§^,^Kk^^M-^t) = ^<ip) E^±/r.±'" • (6-7) 



/ 



One readily recognizes the form of the Kadanoff-Baym ansatz here. Thus, our effective 
Wightman functions can be understood as an en/argerf Kadanoff-Baym ansatz, which allows 
to encode the coherence information as well. Conversely, it is now clear why making the 
KB-ansatz for the propagator from the outset completely misses all coherence information. 

6.2 cQPA Feynman rules 

We will now present the complete set of momentum space Feynman rules to compute the 
effective self-energies iTiQff(k^t) and zneff(A:,t) involving the flavour- coherent propagators 
in the Yukawa theory with interaction (6.1). Let us first write down the general relations 
between the primary CTP-propagators a, 6 = ±, appearing in generic loop-diagrams 
in the CTP-formalism, and the pole- and the Wightman-propagators: 

Q++ = Qt ^QT ^Q< ^ Q+- ^ :^Q< ^ 

G— = = -G^ T G< , G-+ = G> , (6.8) 

where the upper — sign refers to fermionic correlator and lower + sign to scalar corre- 
lator The cQPA Feynman rules can now be stated as follows: 



Draw all perturbative self-energy diagram(s) according to standard CTP-formalism, 
associate the overall signs and symmetry factors for the diagrams and the CTP-indices 
for the propagators. 

For each vertex, except the out-vertex of zSefr(A:,t) or znefr(A:,t), associate a normal 
vertex Feynman rule. For example, in the Yukawa theory (6.1) the vertex rule for 
the incoming scalar line is (here a is the CTP branch index of the vertex): 

yijcPi^PiiPj : -i yij {2n)^S\k' -k+p)a. (6.9) 

The rule for the scalar out- vertex is identical to Eq. (6.9) with complex conjugated 
coupling constant ^ ^ y*. The out- vertex of the self-energy has the same rule but 
without 4-momentum conservation delta function. 



Use the relations (6.8) to write the CTP-propagators in terms of zG^'" and iG 

For all pole-type propagators iG'^'^'*°'*o(fc) substitute the standard (vacuum) propa- 
gator {e.g. Eq. (2.36)) and integrate over k : J ^^^^4 as usual. (Here iG^^'^^ refer to 
the pole-parts of the full Feynman and anti-Feynman propagators.)^^ 



^°Sometimes it may be useful to include the vacuum parts of the Wightman functions iG^'^ to the 
Feynman and anti-Feynman propagators, and use iG^^'^^ instead of zG^'" as the basic pole-propagators. 
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Figure 3: The Feynman rules for the pole-type propagators iQ^<<^'to,to g^j^^ tj^g scalar in- vertices. 



a, I 



k' 



'0' 



•^ti{k) F^if{k, k'; t) A%{k''^ 



a, k 



b,l 



Figure 4: The Feynman rules for the composite Wightman propagators iG^'^ . 

• For each Wightman function substitute a composite propagator iG^^ {k^k') given 
by Eqs. (6.4) and integrate over both 4-momentum variables k and k'\ J ^^^^4 ^^^^4 . 

These rules are presented graphically in the figures 3 and 4, where the arrows in the 
propagator lines indicate the direction of the flow corresponding to the 4-momentum of the 
positive energy state. Although these rules are specifically designed to compute the self- 
energy functions iTiQf^{k,t) and ineff(/c,t), there is no restriction for computing arbitrary 
diagrams with arbitrary number of external legs however. In this case, one would need to 
extend the rules to include in- and out-propagators as well. More details and examples 
related to cQPA Feynman rules in flavour- covariant formulation are presented in [21]. 

7. Application to CP- violating flavour mixing in the presence of collisions 
7.1 Setup 

As an application to our formalism, we consider a scenario motivated by electroweak baryo- 
genesis (EWBG): N x N fermionic mixing through a slowly- varying time-dependent and 
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CP- violating mass matrix, in the presence of Yukawa- interactions with thermal bath. To 
be specific, the interaction Lagrangian is given by 



Ant = -Vi i^ihi^qR + h.C. , (7.1) 

where ijji^ I — 1, . . . , A^, are the mixing fermionic fields in the flavour basis, ifjq and (j) are a 
fermionic field and a complex scalar field assumed to be in thermal equilibrium, and are 
flavour-sensitive Yukawa couplings. First, we observe that the leading order (one-loop) self- 
energies for do not involve out-of-equilibrium propagators and coherence contributions, 
and therefore the effective self-energy functions in Eq. (4.16) reduce to standard Wigner 
representation self-energies i;^'^(/c,t) (see discussion in section 4). Apart from the flavour- 
coupling matrix, these thermal one-loop self-energies are completely identical to the ones 
for single-flavour case, considered in ref. [2]: 



(5]<7'-S3<k-7)i^L, (7.2) 



where iT^Q^{kQ, k, t) are real functions, which eventually need to be evaluated on the mass- 
shells of the mass-basis quasiparticles (see Eq. (5.19) and Eq. (7.10) below): 

iS<(±cj„k,t) = ^-^|/i(±cj„k,t)|, (7.3) 

il^<i±iOi,-k,t) = ±^ i^\hi±ui,-k,t)\ - i±^\Ioi±coi,-k,t)\j , (7.4) 



with 



and 



_ m? + - ml kp _ 

" 2m2 T~ rr^T ' 

_ AV2(mf,m^,m|) |k| _ Vip,|k| 

2mj T ~ mjT ' ^ ' 

where and are the energy and the momentum of the decay products in the decay 
frame, ^/s is the invariant mass of the decaying (heaviest) particle and A(a, c) = {a + 
h — cf' — Abe is the usual kinematic phase-space function. -'^-'^ Furthermore, the thermal 
self-energies S<'>(/c,t) are related by KMS relation [30]: S>(A:,t) = e^^oS<(A:, t), which 
follows directly from the KMS relations of the thermal correlators S^eq and A^'^. 
After rotation to the mass basis, the self-energy is given by 

5]< ik,t) = XS<y t = yl (e< 7° - S< k • 7) Pl , (7.7) 



^The time-dependence in Eqs. (7.4-7.5) comes through time-dependent masses rrii — mi(t). 
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where the rotated couphng matrix is yf- = {y^ij = Unyiy'lfUl,-. Using Eq. (7.7) together 
with the KMS relation for the effective self-energies iiy^^{k^t) it is easy to compute the 
flavour structure of collision terms from (4.21 - 4.22) to get: 

' ([^SL, ^(1 T hv^) f-<]- + [iS>^, 1(1 ± hv^) (7.8) 



mtj 



T ^ ^(1 T hv,,) /-];. + [is>^, 1(1 ± hv^) /-<];.) (7.9) 



1 /Ac^ , /i|k| \ Uij (\-y> m fc<^c r.^> rn nm<^c \ r . . 



where the helicity-projected scalar self-energy functions are given by 

zE>.^±(k,t) = ^(zE>(±a;,-,k,i) + /iiS3>(±a;„k,i)) • (7.10) 

To proceed further we need to specify the flavour-mixing scenario. We shall restrict 
ourselves to the case of two flavours. The schematic structure of the dispersion relations 
for two-flavour mixing is presented in Fig 2. We use the following parametrization for the 
time-dependent flavour non-diagonal mass matrix: 



rriA v(t)e-'<^^ 

^(t)e^^W TUB 



where the off-diagonal magnitude v{t) and the phase a{t) are assumed to have the following 
functional forms motivated by EWBG considerations (see ref. [14]): 



„,)^|(l.tanh(±)). 

.(,) = ^(l + ta„h(A)). (r.l2) 



We have chosen here a hermitian mass matrix for simplicity, such that it is diagonalized 
by a unitary transformation: 

= (/mt/t Uit) = ( ~:f<'' \ . (7.13) 

V m2 / ^' I -sin6'(i)e''"(*) cos6'(i) / 



with 



tan(2^(t)) = ^""^^^ . (7.14) 
rriA - ms 
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Using Eq. (7.13) with y = (7 it is now easy to compute the mixing-gradient matrices (3.20) 
appearing in the Boltzmann equations (3.14-3.15): 

while = 0. We also notice that initially flavour and mass basis coincide. 

Next, we solve the Boltzmann equations (3.14 - 3.15) numerically, using the self- 
energies (7.10) in the collision integrals (7.8 - 7.9) and Eq. (7.15) for the mixing-gradient 
matrix. We assume that the distribution functions fa are initially at tin = — 50/r in 
thermal equilibrium with vanishing chemical potential: 

/|-^±(k, tin) = ±!^fe^{±UJi)dij , f^j^^iK tin) = . (7.16) 

If not stated otherwise, we take = 10/T for a representative transition time-scale. For 
the mixing parameters we use vq = T, ctq = 7r/2, rriA = 2.2 T, uib = 2.0 T, and for the 
interaction parameters we choose massless thermal background fields: rriq = m^^ = 0, and 
the flavour-basis coupling constants y^, I — A^B^ are chosen to be = 1 and yu — 0.5. 
Note that the quasiparticle excitations have time-varying interaction strengths with the 
thermal background: yfiit) and ^22(^)5 which are in general different. 

We solved the system for 60 momentum |k| bins between 10~^T and 30 T, with larger 
bin density in the IR-region of the phase space. The full fermionic correlator has 8 real 
diagonal degrees of freedom and 12 complex non-diagonal degrees of freedom, which each 
has an independent phase and amplitude. The overall discretized numerical problem thus 
consists of 1920 real coupled ordinary differential equations. 

7.2 Numerical results 

Let us now study different observables resulting from this sample computation. In flgure 5 
we show the evolution of the total excesses (over thermal equilibrium) of combined particle 
and antiparticle number densities of mass-states i = 1,2, given by: 

5A^^(|k|,t) = ^ ^(n^k/i + ^zk/i)(^) - (^ik + ^ik)eq , (7.17) 
h 

where the number densities n^k/i and n^k/i are related to flavour-diagonal distribution func- 
tions by Eq. (2.26). Left panel shows the excess 6N1 in heavier fermion species 1 and the 
right panel the corresponding deflciency —SN2 of the total particle number for the lighter 
species 2. The particle number is created at the transition time as a result of coherent 
mixing and interactions, and it is strongest at small momentum region |k| ^ T. The heavy 
particle number is depleted by decays when tT > 200. The deflciency in the lighter particle 
sector is depleted by inverse decays about four times slower, because the particle 2 has a 
weaker Yukawa coupling to the thermal bath. 

A similar model for scalar fields has been considered recently in ref. [14]. However, 
the formalism used in [14] only accounts for the direct (particle-particle or antiparticle- 
antiparticle) flavour-mixing but neglects all particle-antiparticle correlations. The authors 
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Figure 5: Shown is the excess total particle-number density 5Ni for the heavy field i = 1 (left) 
and deficiency of the total particle- number density —5N2 for the light field i = 2 (right). 

of ref. [14] argued that the latter are irrelevant in slowly-varying backgrounds, as their 
effects should average out because of the high-frequency (^ 2cj) oscillations in this sector. 
Since our formalism includes particle-antiparticle coherence (/'^-solutions), we can study 
this phenomenon in detail. To this end, we performed a series of computations both 
including and neglecting the /'^-solutions and compared the results. Somewhat surprisingly, 
we find that the /'^-correlations cannot in general be neglected even in the case of relatively 
slowly varying mass profile with 2^ (10 — 20)/T, and moreover, the particle production 
induced by f^-teims can be the dominant contribution to CP- violation in spontaneous 
baryogenesis models. 

The reason behind this observation is that the relevant scale for the /'^-induced par- 
ticle production is the gradient scale ^ and not the zitterbewegung scale ^ 2lj as 
one might naively expect. Indeed, the structure of the Boltzmann equations (3.14 - 3.15) 
with several cross-couplings between the on-shell functions fa implies that the /'^-solutions 
roughly consists of two qualitatively different modes: the homogeneous mode oscillating 
at the frequency ^ 2a;, and the non-oscillatory forced mode due to the "inhomogeneous" 
(gradient) cross-coupling terms to thermal bulk of the particle excitations and to pertur- 
bative vacuum. While it is true that the homogeneous oscillating modes of /'^-solutions 
are barely excited at all, non-negligible forced-mode excitations are typically generated 
in the transition region with non-vanishing mass- gradients, which give rise to substantial 
particle-antiparticle coherence. 

To display the effect of /'^-correlators we first consider the solutions for the particle 
and antiparticle flavour- coherence correlators |/]^^^(k, t)|. These solutions have nontrivial 
characteristics in the presence of mixing and CP- violation, when the time evolution of the 
system is computed both with and without the particle-antiparticle /'^-correlators. In Fig. 6 
we show the evolution of this correlator for helicity /i = +1 for the particles (left panels) 
and for antiparticles (right panels). Upper panels show the results of a full calculation using 
our complete quantum Boltzmann equations and lower panels the results of a calculation 

^^Here homogeneous and inhomogeneous terms or solutions refer to classification in differential equation 
theory, and not to the spacetime symmetry properties of the system. 
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Figure 6: Shown are the amphtudes of particle (left panels) and antiparticle (right panels) flavour- 
coherence correlators |/]^^_^(k, t) | for positive helicity h = 1. The upper panels are computed with 
the full quantum Boltzmann equations, while the lower panels are computed with the restricted 
equations where the particle-antiparticle /^-correlators are neglected throughout the calculation. 




Figure 7: Shown are the real (left panel) and imaginary (right panel) parts of the correlator 
/i2/i+(^'^) positive helicity = 1, corresponding to the upper left panel in Fig. 6. 



where the /'^-correlators have been forced to zero, emulating the calculation performed in 
ref. [14]. (Note that the amplitude in this figure, despite the same color coding, is two orders 
of magnitude smaller than in Fig. 5.) At small momenta |k| < T, the results of the full 
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Figure 8: Shown are the ampUtudes of the particle- antiparticle flavour- coherence correlators for 
positive helicity h = 1. The flavour-diagonal components |/fi^±(k, t) \ for the heavy field i = 1 (upper 
left) and the light field i = 2 (upper right), and the flavour off-diagonal components \f^2h±0^^^) \ 
state (+) (lower left) and state (— ) (lower right). The results are computed with the full quantum 
Boltzmann equations. 

and restricted calculation are very similar; the normal flavour-mixing dynamics dominates 
the correlator at small momenta. However, in full solution we see also a very wide band 
of excitations at |k| ^ T, which are completely missed in pure flavour-mixing dynamics 
without particle-antiparticle coherence. If one takes the difference of the two results, one 
further finds that the /'^-induced particle-production effect is essentially restricted to large 
momenta. At small momenta |k| < T the residual effect has about tenth of the amplitude 
compared to the large momentum region. Indeed, it is clear that the importance of this 
contribution may only come from its extension over a very large region in phase space. 

So far we have shown only amplitudes of the complex off-diagonal correlators. In these 
plots the oscillatory behaviour of these functions is mostly averaged out. In order to better 
show the true oscillatory structure, we display in Fig. 7 the real and imaginary parts of the 
correlator /i2/z+(k,^), corresponding to the case shown in the upper left panel in Fig. 6. 
The small-scale oscillations, clearly visible both in time and in momentum variable, make 
this type of calculations very challenging numerically. 

In figure 8 we show the particle-antiparticle /'^-correlators for the full calculation. 
The upper panels show the flavour-diagonal correlators |/^^^^^(k, t)| with the heavy flavour 
z = 1 on the left and the light flavour z = 2 on the right, and the lower panels show 
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Figure 9: Shown is the amphtude of particle flavour-coherence correlator |/]^^^(k, t)| for positive 
helicity = 1 for a faster transition with = 5/T (left) and for a slower transition with = 20/T 
(right). The results are computed with the full quantum Boltzmann equations. 



the flavour off-diagonal correlators \ f^2h±(^^^)\- Note that the latter represent coherence 
between particle and antiparticle states of different flavours, i.e. the particle-antiparticle 
flavour coherence. As one would expect these correlators are rather well confined to the 
transition region. However, in this region they are equally large as the flavour off-diagonal 
correlators in particle and antiparticle sectors for large momenta, and they extend much 
farther in momenta than any of the correlators in the restricted calculation neglecting f^- 
solutions. Clearly, these particle-antiparticle correlators, albeit of their short duration, act 
as catalysts at the transition time, giving rise to the large-momentum contribution to the 
flavour off-diagonal /^-correlators. 

Taking a closer look at the full solutions in Fig. 6, we observe that the large momentum 
excitations peak at around |k| ^ 15T. The origin of this enhancement is a result of a 
delicate chain of interactions in the full equation network. However, one can see that this 
region roughly corresponds to momenta for which the off-diagonal particle or antiparticle 
flavour-oscillation frequency equals to the inverse transition width^^: ^ Acjij. In this 
case the flavour oscillations are resonant with the rate of change induced by the background 
variation in the transport equation (3.14), which can lead to large enhancement of the 
flavour off-diagonal correlator /j^/^^. Note that this enhancement works typically only 
for the UV- modes, since the resonance condition yields UV-momenta for relatively large 
transition widths r^. With the parameters used in the current example we have 



|k| 



9 9 

nil ~ ^2 
2T2 



T^T ^ 4-40, (7.18) 



where the given range results from the variation of the mass difference m^(t) — m2{t) over 
the transition region as v{t)/T evolves from zero to one. 



Remember that while particle and antiparticle flavour-oscillation shells correspond to ko = icjij the 
characteristic oscillation frequencies of these correlators are :h2Aujij. The situation is inverted for the 
particle-antiparticle flavour coherence, living on the shells ztAcj^j but oscillating with frequencies itcjij. 
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Figure 10: Shown are the charge asymmetries j^(|k|,t) for heavy field i = 1 (left panels) and for 
light field i = 2 (right panels). The upper panels are computed with the full quantum Boltzmann 
equations, while the lower panels are computed with the restricted equations where the particle- 
antiparticle /'^-correlators are neglected throughout the calculation. 



To test the above explanation we have redone our full calculation with two other 
transition times and show the corresponding flavour off-diagonal correlators in Fig. 9. Left 
panel corresponds to a relatively rapid transition with = 5/T and the right panel to a 
much broader one with = 20/T. In the former case the resonant enhancement region is 
narrower, as expected on the basis of Eq. (7.18). In the latter case the overall amplitude is 
getting very weak. However even in the latter case the result differs significantly from the 
one where the particle-antiparticle correlators are neglected. 

Finally, let us consider the quantity which is most interesting for many applications, 
the produced charge asymmetry: 

(|k|,t) = 5^(n,k/z - n,k/^)(t) . (7.19) 

h 

We plot this quantity in Fig. 10, again for the heavy fields (left panels) and the light 
fields (right panels) and with full (upper panels) and restricted (lower panels) equations of 
motion. Even for fairly broad transition considered here (we again use = 10/T), there 
is a considerable asymmetry for near UV modes above k > 3T which is completely missed 
in the restricted calculation where /'^-correlators are neglected. Because of the large phase 
space, this effect dominates the total asymmetry generation. To make this clearly visible. 
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we consider the time evolution of the total asymmetry integrated over the phase space: 



fit) = = / «^J?(|k|,i) . (7.20) 

i i 

We show this quantity in Fig. 11 for three different transition times: = 5/T, lO/T 
and 20/r. The solid curves correspond to the full calculation and the dashed ones to 
the restricted calculation. In particular for small transition times the influence of the 
particle-antiparticle correlators is crucial and the effect dominates over the pure flavour 
oscillation mechanism. Even for rather smooth transitions = 20/T (lowest curves) the 
effect is quantitatively significant. We hence conclude that neglecting particle-antiparticle 
correlators is in general not warranted. The late time tT > 100 behaviour of the asymmetry 
is highly dependent on the type of interactions in the model. Using only three-body 
(1^2) interactions in our toy model leads to a significantly suppressed UV-mode charge 
transfer rate with the external field tpq^ especially for the light field 2, leading a to maximal 
asymmetry at very late times tT ^ 400 — 600 and to a very slow chemical equilibration 
rate, for which an exact result can be derived from equations (7.4 - 7.10): 

rchem,(|k|) = ^^--ln( ^ r,,em,2||kM0T^4xl0-4T. (7.21) 

6Z7r LOi |K| sinh( ) 

This rate agrees with the ordinary one particle decay rate with the expected order of 
magnitude rchem,2 ~ 10~^T in the thermal region |k| ~ T. The very slow integrated 
rate seen here reflects the fact that the late time decays are dominated by the weakly 
interacting light state 2 and that the asymmetry is mostly confined to UV modes at |k| ^ 
15 — 20 T. Taking into consideration also the four-body (2 ^ 2) interactions should make 
the thermalization process much faster especially for UV modes. However, the considerable 
difference in asymmetry generation in the transition region tT < 50 between the full and 
the restricted cQPA calculation persists regardless of the model-specific interaction details. 

Let us stress that the /'^-induced enhancement of the asymmetry generation is an 
entirely new type of phenomenon, which would not have been discovered without full 
cQPA transport formalism. We expect this effect to be important for many applications, 
and in particular to the electroweak baryogenesis type scenarios. 

8. Flavour coherence from the operator formahsm 

In this section we use standard operator formalism of QFT to study flavour-mixing com- 
plex scalar fields in constant background. We show that flavour coherence is emergent 
provided that the off-diagonal (i 7^ j) expectation values {oli^cijk')^ (^Ik^ikO or (^zk^k') 
are nonvanishing. In this case, and further satisfying spatial homogeneity and isotropy, we 

^^More precisely, this quantity corresponds to the expectation value of the number current density: 
— ^i{'4^il^^i) — f (f:;^Tr[5'^(/c, t)], where the vacuum part is subtracted and the trace is over 
both Dirac and flavour indices. Note that this quantity is flavour-basis invariant. 
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Figure 11: Shown is the evolution of the total asymmetry integrated over the phase space j^{t) = 
Ji (^)+J2 (^)- The solid lines correspond to the solution of the full equations of motion and the dashed 
lines to the restricted equations of motion with = 0, both respectively with rw = 5/T, 10/T and 
20/T from top to bottom. 



then find that the corresponding 2-point correlator A^(/c,t) has the cQPA-shell structure 
with ko = zbo; or ko = zbAcj. We show that the desired nonstandard expectation values 
can be realized by certain class of exponentiated superposition states, flavoured squeezed 
states, obtained from the vacuum state by unitary squeezing operation [24]. We also show 
that flavoured squeezing corresponds to flavoured Bogolyubov transformations, providing a 
generalization of the standard non-flavoured case. We further present a one-to-one identifi- 
cation between the on-shell functions f^J^ and the independent parameters of the flavoured 
squeezed state or the corresponding flavoured Bogolyubov transformation. 

We consider mixing complex scalar fields 0^, i = 1,...,A^, with standard operator 
expansions 

where cj^k = y^k^ + mf and a^^ (^ik) and b^^ (6^k) are the creation (annihilation) oper- 
ators for particles and antiparticles of mass eigenstates i = 1, . . . , A^, respectively. These 
operators satisfy canonical commutation relations 

[aik,a^v] =[4k'4k^] = [^ik,&jk^] = [^^k'^^k^] = 0' 

[«.k, =[&.k, &]k^] = 2c^,k(2^)'5'(k - k')S,j , (8.2) 

and all commutators between a^k and bj\^f (and their hermit ian conjugates) are vanishing. 
Let us now assume a state with the expectation values 

= + /^_<(|k|)) 20ijU27rf6'{k - k') , 
{hi. a,v) = (a]k, 4)* = /;=<+(|k|)) 2cD,,k(2vr)353(k + k') , (8.3) 
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consistent with the requirement of spatial homogeneity and isotropy. Using the operator 
expansion (8.1) and the relations (8.3) it is straightforward to compute the Wightman 
function: 



i^u^v) ={4>\{v)ct>i{u)) 



-I 



d^k 



ij ^zk-(u-v) 



fm< -i{uJiUo—ujjVo) _ fm< i{uJiUo-ujjVo) 
Jij-\- ^ Jij- ^ 

\fc< -i{uJiUo+ujjVo) I fc<* i{uJiUo-\-ujjVo) 



(8.4) 



where we have dropped the k-subscripts of cjjk and Wjjk- The Wigner transformation then 
gives {t = XQ^ (no + vo)/2) 

iAfj{k,x) = j d^(u-v)e^'=o("°-^°)-*'^-("-^)iA,^.(u,?;) 

+^f+e-*2-*5(fco - Aa;) + fJ<*e'^'''S{ko + Au)] , (8.5) 
which is readily identified as the singular cQPA Wightman function of Eq. (5.6) by 

/-t<(k, t) = ^7<(k) e^^^A.,, t ^ fc<^^^^ ^) ^ ^< (k) e^^^c.., t ^ (8.6) 



fij- 



using the relation = fj^- Indeed, we have found that the correlator A<(/c, t) involves 
flavour coherence and all four singular shells ko = ±0 and ko = dzAo; are excited, provided 
that the expectation values of Eq. (8.3) are nonvanishing (and consistent with the spatial 
homogeneity and isotropy). 

Next, we would like to find a state with the nonstandard expectation values (8.3). 
It is clear that any number state: (H ^J^)(n with a finite number of quanta cannot 

m n 

possibly satisfy the relations (8.3), because the flavour off-diagonal expectation values 
would necessarily vanish. It turns out, however, that certain superposition states with 
unspecified (infinite) number of quanta indeed have the desired properties (8.3). To see 
this, we use the methods of squeezed state formalism [24]. 

We begin by defining a unitary two-mode squeeze operator for multiple mixing fields: 



S2{rij) 



— — 



exp 



(8.7) 



where rij(|k|) is a complex matrix in flavour indices, and we denote ajk = a^k/ V'^'^ik ^^'^ 
bii^ = 6ik/\/2tJik- Next, we show how the creation and annihilation operators transform 
upon acting on S2, i.e. we want to compute S'2aik5'2 and its hermitian conjugate. First, 
we find that 



(8.8) 
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Furthermore, r^j can be diagonalized by a biunitary transformation: r — U'^r^V with 
(^d)ij — "Ti^ij^ where n(|k|) > and C/(|k|) and l^(|k|) are unitary matrices. This suggests 
to define the following linear combinations: 



(8.9) 



which by construction satisfy simple diagonal commutation relations 

[Qk±,^] = TnQk±. (8.10) 

By recursive use of these relations we find that Ci\^±A^ = (A =f r^)^Qkdz5 furthermore 

52C,k±5| = e±^^c,k±. (8.11) 

From this equation and the inverse of Eq. (8.9) we directly obtain the transformation law 
for ai\^: 

j 



j 

= ^(Ajaik + 7ij^j-k)' 



.12) 



with 



«i,(|k|)^[Wcosh(rd)F]^. 
7i,(|k|)^[t/tcosh(rd)C/].. 



= [uUmh{ra)V] 



(8.13) 



We see that the transformation law (8.12) has the form of a flavoured Bogolyubov transfor- 
mation. For the inverse transformation we get 



j 



(8.14) 



Because of the unitarity of the transformation (8.12), it is clear that the transformed 
operators a^k? ^ik satisfy the same commutation algebra as a^k? ^ik- Also, we instantly note 
that the state |0) = 5210) is a vacuum state for the transformed operators, i.e. a^k|0) = 0. 

It is now easy to construct a state |Q) with the desired properties (8.3). It turns out 
that a first guess: 

|0) (8.15) 



\n) = |0) = S2\0) = exp l^^l ^ {riMj_^ - rtja]^bl^) 
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is indeed correct, since a direct computation gives 



= .. {27rfS^{k - kO (8.16) 

and similarly 

(4 ^ik') = 2VS^ (27r)3(53(k - k') , 

(&ik a,v) = -(/3^),i 2./m] (27r)3(53(k + k') . (8.17) 

These relations readily provide a one-to-one identification between the on-shell functions 
and the parameters of the flavoured Bogolyubov transformations (or equivalently the pa- 
rameters Tij of the corresponding squeezed state): 

9. Conclusions and outlook 

In this paper we have derived quantum transport equations and perturbative Feynman 
rules for flavour-mixing quantum fields in spatially homogeneous systems, including both 
flavour- and particle-antiparticle coherence. We have considered both fermionic and bosonic 
fields using flavour-mixing Yukawa- couplings as a model for interactions. Our formalism is 
based on the coherent quasiparticle approximation (cQPA) [1-6,21], which is a distributive 
expansion of the 2-point correlation functions in the limit of small background gradients and 
weak interactions. In cQPA the coherence information, both of the flavour- and particle- 
antiparticle mixing, resides on distinct coherence shells in the phase space, located at 
time-like average energies ujij = ^{uji + uj) for direct flavour- mixing and at space-like 
average energies (between positive and negative energy states) Aojij = ^{uji — uoj) for 
particle-antiparticle flavour-mixing. 

We have derived explicit, generic forms of flavour-coherent transport equations for 
fermion and scalar fields, which are valid for any type of interactions. These transport 
equations reduce to the usual quantum Boltzmann equations when the coherence correla- 
tors are neglected. We have also derived a simple set of Feynman rules for perturbative 
calculations of the self-energy functions appearing in the collision integrals of the transport 
equations including non-local coherence effects. For the single-flavour case the Feynman 
rules derived in this paper are equivalent but simpler than the ones introduced in ref. [6]. 
In a companion paper [21] we present an alternative, flavour-covariant formulation of the 
cQPA transport equations and Feynman rules. This formulation also provides a straight- 
forward generalization of the formalism to the case of nonzero dispersive self-energy S/^ 
with nontrivial quasiparticle dispersion relations. 



fC< _ 



I J 1 



Jij- 



UJjUJj 



^.18) 
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As a numerical example, we have studied a two-flavour mixing scenario, where the 
mixing fermions interact with a thermal bath. A similar model for mixing scalar fields was 
recently considered in ref. [14], where the particle-antiparticle coherence solutions were 
neglected in the dynamics based on physical arguments. However, in our calculations 
we have found that neglecting particle-antiparticle coherence is not in general justified 
even when the background field is slowly varying in comparison with the zitterbewegung 
frequency ^ 2bj. Indeed, particle-antiparticle coherence correlators were explicitly shown 
to be excited in such cases. Moreover, these correlators can trigger a resonant growth 
of flavour coherence at large momenta over a large part of the phase space. Comparing 
the full calculations with the restricted model calculations, where the particle-antiparticle 
coherence correlators were neglected, we have found that these new modes can dominate 
over the pure flavour dynamics. These results may suggest a new, more efficient way of 
generating an asymmetry in the EWBG-type models, even for relatively slowly-varying 
mass profiles such as in the case of electroweak phase transition in MSSM. 

We have also shown that cQPA coherence correlators are related to squeezed states 
in the operator formalism language. Indeed, we have found that the on-shell coherence 
functions cannot be related to any number states, but a one-to-one connection exists be- 
tween the coherence functions and the parameters of a squeezed state. We also showed that 
flavoured squeezing corresponds to flavoured Bogolyubov transformations, generalizing the 
usual squeezing formalism to the flavour-mixing case. 

Our formalism has several interesting applications. First, a study of thermal leptogen- 
esis in the resonant regime using cQPA formalism is in progress. Second, we are working on 
to derive transport equations for the usual neutrino-mixing scenarios from first principles 
using our formalism. Third, it will be interesting to generalize the simple toy-model for 
the quenching EWBG-transition introduced in [6] to the more realistic flavour-mixing case. 
Fourth, we are working on to extend our results to stationary, planar symmetric problems, 
which are more directly relevant for traditional EWBG models. 
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A. Gradient corrections to the phase space structure of A^'^(/c,t) 

In this section we examine the first order gradient corrections ~ dtm to the phase space 
structure of the scalar Wightman functions A'^'^{k,t). We also find generalizations of the 
relations (5.12) between the moments pn and the on-shell functions fa, which can be used 
to obtain order 0(rdtm) corrections to the collision integrals (5.20-5.22). 

Using both KB-equations (5.2-5.3) and keeping the first order mass gradients we obtain 
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a generalization of the (algebraic) zeroth-order constraint equation (5.4): 



=ko {{kl - - M^)kl + (Ammf) iA^'> 



d 



- [kQ^Auu + u'Auj) - koAuju;{AuAuj + uCj)) dk^A^ 



(A.l) 



Because of the fco-derivatives dkoA^''^ {k,t), this equation cannot be solved by a simple 
singular phase-space structure like Eq. (5.6). To proceed, we integrate over ko to obtain 
the following relations for the moments pn- 

pn = pn-2{AuJ^ + O)^) - pn-^Auj'^U)'^ 

i(2n — 5) / . / / . X i(2n — 7) , , . , ^ / x / * 
+ ^ — ^ -pn-siAuj'Cj + cj'Acj) - ^ — -pn-5AujCu{Auj'Auj + Cu'o) . (A.2) 

Clearly, the moments of the singular (zeroth order) correlator (5.6) satisfy the zeroth-order 
terms of this equation. To solve Eq. (A.2) up to first order, we make an ansatz for the 
Wightman functions A^''^, where a linear correction is parametrized as an expansion in 
derivatives of (zeroth order) phase-space delta functions: 

A<'> (A;, t) = A<^ {k,t) + J2 91 {±an±5{ko T ^) + bn±S{ko T Ac^)) , (A.3) 

n± 



where A^q {k,t) is the zeroth-order solution (5.6), and and bn± are of order dtm. 
Parameters an± and bn± can be determined by computing the moments pn from the ansatz 
(A.3) and using them in Eq. (A.2). We find that only first- and second-order /cQ-derivatives 
of the delta functions are required to parametrize the first-order gradient corrections: 



UJjUJ 



J ± L 



i Auj' uj — uj' Auj ^ , i ^ . . 



c<,> 

± 



i Auj'uj — uj' Auj 



(A.4) 



As a distribution, i.e. if Eq. (A.4) is multiplied by a smooth test function, it is possible 
to replace the derivatives of the delta-functions by a Gaussian distribution function with a 
complex width and a shift in the expectation value. This can be shown straightforwardly 
by comparing the resulting moment-integrals up to the first-order terms. Therefore, we 
can parametrize the first-order Wightman function equivalently as 



E 



'J ± 



I nm<,> 



, _ i Auo'uo — uo' Auj i . , 



fco K r f . X i Auj'u — u' Au , i , 



where 



exp 



2 cD2 - Aa;2 
2a2 



1-j 



(A.5) 



(A.6) 
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is the (Gaussian) normal distribution function. According to the expectations, we have 
found that the first-order gradient terms give rise to a finite width for the Wightman 
functions A^'^(A;,t), which is directly proportional to dtm. Also, a shift to the positions 
of the "shells" (also proportional dtm) is induced. 

The extended phase-space structure (A. 4) (or (A. 5)) for A^'^(fc,t) allows us to com- 
pute gradient corrections to the relations (5.12) and consequently to the collision integrals 
(5.20-5.22). First, by using Eq. (A. 4) in the kinetic equation (5.3) and taking any four 
independent moments (or equivalently using the moment equations (5.14)) we obtain the 
equations 



m<,> 



/r '> = 0, 



C<,> 



/r'^ = o. 



(A.7) 



cD(a)2 - Aa;2) 

which are valid up to order 0{dtm) (here we neglect the gradients of the mixing matrix U 
for simphcity). By taking the moments (A. 4) and using Eqs. (A.7) to relate dtpo to fa we 
then obtain the following generalization of the invertible relations (5.12): 



/ P. \ 

dopo 
Pi 

\ J 



13 



( 



1 

-2iAa; 



-1 
-2iAa; 



1 1 

— 2iw liuj 
Au) -Alo 
Au^ Au^ 



+ X 



jr. 

n 
\f-J 



(A.8) 



IJ 



with 



Xij — 



( 

\ + \AJ 





2L>2W + \AJ 





2D2AUJ + 



where 



(jo'uJ 



Aoj'Au 



and 





-2D2ALO - \uj' 
. Aoj'uj — co'Alo 



(A.9) 



These relations can be used to obtain order OiVdtm) corrections to the collision integrals 
(5.20-5.22). However, the corrections of this order involve also direct derivative terms 
^ dtV Tdtm and the mixing gradient terms ^ = iUdfU'^ ^ as well as corrections to the 
resummation of the oscillatory coherence terms in the collision integrals, which have not 
been considered here. Note that the flow terms of Eqs. (5.14) obtain no gradient corrections 
as they are written in terms of the moments. 



References 

[1] M. Herranen, K. Kainulainen and P. M. Rahkila, Towards a kinetic theory for fermions with 
quantum coherence, Nucl. Phys. 5 810 (2009) 389 [arXiv:0807.1415 [hep-ph]]. 

[2] M. Herranen, K. Kainulainen and P. M. Rahkila, Quantum kinetic theory for fermions in 
temporally varying backgrounds, JHEP 0809 (2008) 032 [arXiv:0807.1435 [hep-ph]]. 



-46- 



[3] M. Herranen, K. Kainulainen and P. M. Rahkila, Kinetic theory for scalar fields with nonlocal 
quantum coherence, J HEP 020b (2009) 119 [arXiv:0812.4029v2 [hep-ph]]. 

[4] M. Herranen, Quantum kinetic theory with nonlocal coherence, PhD Thesis, University of 
Jyvaskyla Res. Rep. 3/2009 [arXiv:0906.3136 [hep-ph]]. 

[5] M. Herranen, K. Kainulainen and P. M. Rahkila, Coherent quasiparticle approximation cQPA 
and nonlocal coherence, J. Phys. Conf. Ser. 220 (2010) 012007 [arXiv:0912.2490 [hep-ph]]. 

[6] M. Herranen, K. Kainulainen and P. M. Rahkila, Coherent quantum Boltzmann equations 
from cQPA, JHEP 1012 (2010) 072 [arXiv: 1006. 1929 [hep-ph]]. 

[7] R. Barbieri and A. Dolgov, Bounds On Sterile- Neutrinos From Nucleosynthesis^ Phys. Lett. 
B 237 (1990) 440; Neutrino oscillations in the early universe, Nucl. Phys. B 349 (1991) 743; 
K. Kainulainen, Light Singlet Neutrinos And The Primordial Nucleosynthesis, Phys. Lett. B 
244 (1990) 191; 

K. Enqvist, K. Kainulainen and J. Maalampi, Resonant Neutrino Transitions And 
Nucleosynthesis, Phys. Lett. B 249, 531 (1990); 

Refraction And Oscillations Of Neutrinos In The Early Universe, Nucl. Phys. B 349 (1991) 
754; 

K. Enqvist, K. Kainulainen and M. J. Thomson, Stringent cosmological hounds on inert 
neutrino mixing, Nucl. Phys. B 373 (1992) 498; 

D. Boyanovsky, C. M. Ho, Charged lepton mixing and oscillations from neutrino mixing in 
the early universe, Astropart. Phys. 27 (2007) 99-112; 

C. M. Ho, D. Boyanovsky, H. J. de Vega, Neutrino oscillations in the early universe: A Real 
time formulation, Phys. Rev. D72 (2005) 085016. 

[8] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, On The Anomalous Electroweak 
Baryon Number Nonconservation In The Early Universe, Phys. Lett. B 155 (1985) 36; 
A. G. Cohen, D. B. Kaplan and A. E. Nelson, Progress in electroweak baryogenesis, Ann. 
Rev. Nucl. Part. Sci. 43 (1993) 27 [arXiv:hep-ph/9302210]; 

V.A. Rubakov and M.E. Shaposhnikov, Electroweak baryon number non- conservation in the 
early universe and in high-energy collisions, Usp. Fiz. Nauk 166 (1996) 493, Phys. Usp. 39 
(1996) 461 [hep-ph/9603208]. 

[9] M. Joyce, T. Prokopec and N. Turok, Nonlocal electroweak baryogenesis. Part 2: The 

Classical regime, Phys. Rev. D 53, 2958 (1996) [hep-ph/9410282]; Electroweak baryogenesis 
from a classical force, Phys. Rev. Lett. 75, 1695 (1995); [Erratum- ibid. 75, 3375 (1995)] 
[hep-ph/9408339]; Nonlocal electroweak baryogenesis. Part 1: Thin wall regime, Phys. Rev. D 
53, 2930 (1996). 

[10] J. M. Cline, M. Joyce and K. Kainulainen, Supersymmetric electroweak baryogenesis, JHEP 
0007 (2000) 018 [hep-ph/0006119]; Supersymmetric electroweak baryogenesis in the WKB 
approximation, Phys. Lett. B 417 (1998) 79, [Erratum- ibid. B 448 (1999) 321] 
[hep-ph/9708393]; 

J. M. Chne and K. Kainulainen, A new source for electroweak baryogenesis in the MSSM, 
Phys. Rev. Lett. 85 (2000) 5519 [hep-ph/0002272]. 

[11] K. Kainulainen, T. Prokopec, M. G. Schmidt and S. Weinstock, First principle derivation of 
semiclassical force for electroweak baryogenesis, JHEP 0106, 031 (2001) [hep-ph/0105295]; 
Semiclassical force for electroweak baryogenesis: Three-dimensional derivation, 
Phys. Rev. D66 (2002) 043502 [hep-ph/0202177]; Quantum Boltzmann equations for 
electroweak baryogenesis including gauge fields, [hep-ph/0201293]. 



-47- 



[12] T. Prokopec, M. G. Schmidt and S. Weinstock, Transport equations for chiral fermions to 
order h bar and electroweak baryogenesis. Part i., Annals Phys. 314 (2004) 208 
[hep-ph/0312110]; Transport equations for chiral fermions to order h-bar and electroweak 
baryogenesis. Part 11. , Annals Phys. 314 (2004) 267 [hep-ph/0406140]. 

[13] T. Konstandin, T. Prokopec and M. G. Schmidt, Kinetic description of fermion flavor 
mixing and CP-violating sources for baryogenesis, Nucl. Phys. 5 716 (2005) 373 
[arXiv:hep-ph/0410135]; 

T. Konstandin, T. Prokopec, M. G. Schmidt and M. Seco, MSSM electroweak baryogenesis 
and flavour mixing in transport equations, Nucl. Phys. B 738 (2006) 1 
[arXiv:hep-ph/0505103]. 

[14] V. Cirighano, C. Lee, M. J. Ramsey-Musolf and S. Tuhn, Flavored Quantum Boltzmann 
Equations, Phys. Rev. D 81 (2010) 103503 [arXiv:0912.3523 [hep-ph]]. 

[15] V. Cirighano, C. Lee and S. Tuhn, Resonant Flavor Oscillations in Electroweak Baryogenesis, 
arXiv:1106.0747 [hep-ph]. 

[16] B. Garbrecht, T. Prokopec and M. G. Schmidt, Coherent baryogenesis, Phys. Rev. Lett. 92 
(2004) 061303 [arXiv:hep-ph/0304088]; SO(IO) - GUT coherent baryogenesis, Nucl. Phys. B 
736 (2006) 133 [arXiv:hep-ph/0509190]. 

[17] M. Fukugita and T. Yanagida, Baryogenesis Without Grand Unification, Phys. Lett. B 174: 
(1986) 45. 

[18] A. Pilaftsis and T. E. J. Underwood, Resonant leptogenesis, Nucl Phys. B 692 (2004) 303 
[arXiv:hep-ph/0309342]; Electroweak- scale resonant leptogenesis, Phys. Rev. D 72 (2005) 
113001 [arXiv:hep-ph/0506107]. 

[19] A. Abada, S. Davidson, F. X. Josse-Michaux, M. Losada and A. Riotto, Flavour Issues in 
Leptogenesis, JC^P 0604 (2006) 004 [arXiv:hep-ph/0601083]; 

E. Nardi, Y. Nir, E. Roulet and J. Racker, The importance of flavor in leptogenesis, JHEP 
0601 (2006) 164 [arXiv:hep-ph/0601084]; 

M. Garny, A. Hohenegger, A. Kartavtsev and M. Lindner, Systematic approach to 
leptogenesis in nonequilibrium QFT: vertex contribution to the CP-violating parameter, Phys. 
Rev. D 80 (2009) 125027 [arXiv:0909.1559 [hep-ph]]; Systematic approach to leptogenesis in 
nonequilibrium QFT: self- energy contribution to the CP-violating parameter, Phys. Rev. D 81 
(2010) 085027 [arXiv:0911.4122 [hep-ph]]; 

M. Beneke, B. Garbrecht, M. Herranen and P. Schwaher, Finite Number Density Corrections 
to Leptogenesis, Nucl. Phys. B 838 (2010) 1 [arXiv:1002.1326 [hep-ph]]; 

M. Beneke, B. Garbrecht, C. Fidler, M. Herranen and P. Schwaher, Flavoured Leptogenesis in 

the CTP Formalism, Nucl. Phys. B 843 (2011) 177 [arXiv: 1007.4783 [hep-ph]]; 

A. Anisimov, D. Besak and D. Bodeker, Thermal production of relativistic Majorana 

neutrinos: Strong enhancement by multiple soft scattering, J CAP 1103 (2011) 042 

[arXiv: 1012.3784 [hep-ph]]. 

[20] A. Anisimov, W. BuchmuUer, M. Drewes and S. Mendizabal, Leptogenesis from Quantum 
Interference in a Thermal Bath, Phys. Rev. Lett. 104 (2010) 121102 [arXiv: 1001. 3856 
[hep-ph]]; Quantum Leptogenesis I, Annals Phys. 326 (2011) 1998 [arXiv:1012.5821 
[hep-ph]]. 

[21] M. Herranen, K. Kainulainen and P. M. Rahkila, Flavour- coherent propagators and Feynman 
rules: Covariant cQPA formulation, arXiv: 1108.2371 [hep-ph]. 



-48- 



[22] J. F. Koksma, T. Prokopec and M. G. Schmidt, Decoherence in an Interacting Quantum 

Field Theory: The Vacuum Case, Phys. Rev. D 81 (2010) 065030 [arXiv:0910.5733 [hep-th]]; 
Entropy and Correlators in Quantum Field Theory, Annals Phys. 325 (2010) 1277 
[arXiv: 1002.0749 [hep-th]]; Decoherence in an Interacting Quantum Field Theory: Thermal 
Case, Phys. Rev. D 83 (2011) 085011 [arXiv:1102.4713 [hep-th]]. 

[23] A. Giraud and J. Serreau, Decoherence and thermalization of a pure quantum state in 
quantum field theory, Phys. Rev. Lett. 104 (2010) 230405 [arXiv:0910.2570 [hep-ph]]; 
F. Gautier and J. Serreau, Quantum properties of a non- Gaussian state in the large- N 
approximation, Phys. Rev. D 83 (2011) 125004 [arXiv: 1104.0421 [quant-ph]]. 

[24] C. M. Caves and B. L. Schumaker, New formalism for two-photon quantum optics. 1. 
Quadrature phases and squeezed states, Phys. Rev. A 31 (1985) 3068; 
B. L. Schumaker, Quantum mechanical pure states with Gaussian wave functions, Phys. 
Rep. 135 (1986) 317; 

L. P. Grishchuk and Yu. V. Sidorov, Squeezed quantum states of relic gravitons and 
primordial density fluctuations, Phys. Rev. D 42 (1990) 3413; 

A. Albrecht, P. Ferreira, M. Joyce and T. Prokopec, Inflation and squeezed quantum states, 
Phys. Rev. D 50 (1994) 4807 [arXiv:astro-ph/9303001]. 

[25] J. S. Schwinger, Brownian motion of a quantum oscillator, J. Math. Phys. 2 (1961) 407; 
L. V. Keldysh, Diagram technique for nonequilibrium processes, Zh. Eksp. Teor. Fiz. 47 
(1964) 1515 [Sov. Phys. JETP 20 (1965) 1018]. 

[26] E. Calzetta and B. L. Hu, Nonequilibrium Quantum Field Theory, Cambridge, UK: Univ. Pr. 
(2008) 535 p. 

[27] L. Kadanoff and G. Baym, Quantum Statistical Mechanics, Benjamin, New York (1962). 

[28] J. M. Luttinger and J. C. Ward, Ground state energy of a many fermion system. 2, Phys. 
Rev. 118 (1960) 1417; 

J. M. Cornwah, R. Jackiw and E. Tombouhs, Effective Action For Composite Operators, 
Phys. Rev. D 10 (1974) 2428; 

K. c. Chou, Z. b. Su, B. 1. Hao and L. Yu, Equilibrium And Nonequilibrium Formalisms Made 
Unified, Phys. Rept. 118 (1985) 1; 

E. Calzetta and B. L. Hu, Nonequilibrium Quantum Fields: Closed Time Path Effective 
Action, Wigner Function and Boltzmann Equation, Phys. Rev. D 37 (1988) 2878 . 

[29] P. Danielewicz, Quantum Theory Of Nonequilibrium Processes, i. Annals Phys. 152 (1984) 
239; 

S. Mrowczynski and U. W. Heinz, Towards a relativistic transport theory of nuclear matter. 
Annals Phys. 229 (1994) 1 . 

[30] R. Kubo, Statistical Mechanical Theory Of Irreversible Processes. I, J. Phys. Soc. Jap. 12 
(1957) 570; 

R. Kubo, M. Yokota and S. Nakajima, Statistical Mechanical Theory Of Irreversible 
Processes. II, J. Phys. Soc. Jap. 12 (1957) 1203; 

P. C. Martin and J. S. Schwinger, Theory of many particle systems. I, Phys. Rev. 115 (1959) 
1342. 



-49- 



